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PREFACE 


The Electric Punched Card Method, known also as the 
Hollerith, or Electric Accounting Machine Method, is now 
generally known to the business and statistical worlds. Its 
use in physical science has, however, been confined largely to 
astronomy. 

The purpose of this book is to show the px ssibilities of the 
method in scientific computation. It is not a manual of 
operation for the various machines but is designed to enable 
4 scientist so to formulate his problem that any skilled opera- 
tor of the machines could carry it out. 

No account is given of punched card tabulating machines 
in which the cards operate the machines mechanically in- 
stead of electrically. The great advantage of the electric 
method for scientific computation is its flexibility. In the 
electric method the necessary machine changes in going from 
one problem to another can be accomplished by the operator 
in a few minutes. 

Nearly all astr« nomical calculations now being done by the 
Punched Card Method are being performed in The Thomas J. 
Watson Astronomical Computing Bureau, a scientific non- 
profit-making institution which was organized for the use of 
astronomers. Most of the methods in the book will therefore 
be described in terms of the equipment of the Bureau. This 
procedure will assist astronomers in using the Bureau, and will 
furnish for others a working example. A brief description of 
the organization of the Bureau is given in appendix c 

The computing laboratory which is now operated by the 
Bureau was developed by the writer, first in the Columbia 
University Statistical Bureau and later in the Department of 
Astronomy of Columbia University. The directors of these 
institutions, Dr. Ben Wood and Dr. Jan Schilt, have given 
every encouragement. 

The laboratory would not have been possible without the 
interest and generosity of Dr. Thomas fs Watson, President 
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of the International Business Machines Corporation. The 
undertaking has had his cooperation and encouragement from 
the very beginning. 

It is a pleasure to acknowledge the technical and other 
assistance of many officials of the International Business Ma- 
chines Corporation. Mr. G. W. Baehne made every effort to 
secure equipment that would effectively meet all of the re- 
quirements. Mr. William Fitzpatrick, Customers’ Service 
Supervisor of the New York District, assisted by Mr. W. L. 
Seasholtz of the Customers’ Service Division, arranged the 
special circuits for the process of numerical integration de- 
scribed in chapter VIII. The Bureau itself has grown out of 
a suggestion by Mr. C. H. Tomkinson, then Manager of the 
Commercial Research Department. 

Miss Lillian Feinstein has contributed largely to the suc- 
cess of several of the applications in part III. In connection 
with the work on the Yale Zone Catalogues and the Columbia 
Photometric Program she has made valuable suggestions as to 
the arrangement of the work and has personally performed 
most of the computations. 

The author has been able to incorporate in the manuscript 
a number of corrections and improvements suggested by Dr. 
Frank Schlesinger, Dr. Dirk Brouwer, Mr. Hans Hertz, Mr. 
George Essenwine, Miss Lillian Feinstein, Mrs. Grace O. 
Savage and Dr. Robert Wolf. 

Rutherfurd Observatory W. J: sigma 
Columbia University — 


New York City 
December 15, 1939. 
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PART I. THE PUNCHED CARD METHOD 
CHAPTER | 
INTRODUCTION 


There have been three major advances in computing 
methods brought about respectively by the introduction ol 
logarithms, of manually operated computing machines, and 
of automatic tabulators using pun¢ hed cards. The first table 
of logarithms was published by John Napier in 1614, and for 
nearly three centuries all extensive calculations were per- 
formed by this method: The present generation has seen the 
logarithm replaced by hand and electric computing machines. 
These are all similar in that the simple operations of arith- 
metic are performed more or less automatically upon numbers 
which are put into the machine by hand. For the automatit 
tabulators the initial numbers are punched on cards and 
checked. From then on they receive no personal attention, 
as they are automatically transferred thr« ugh flexible wiring 
to the proper place in the machine for the various operations. 

The complete line of equipment using the punched card 
system is too extensive to describe here; suffice it to say that 
there are standard machines for practically all purposes. 
They have been developed for commercial uses but with 
slight modifications become ideal for scientific work. The 
main question in any case is not ‘‘can the problem be s ilved 
by these machines?” but rather “have I enough operations 
of this type or that to justify such powerful equipment?” 
It might be mentioned also that, with standard parts, special 
machines for almost any imaginable purpose could be con- 
structed, but except in very rare Cases the expense and delay 
involved make the standard machines preferable. 

The first extensive use of the early Hollerith Tabulator in 


astronomy was made by Comrie. He used it for building a 
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table from successive differences, and for adding large numbers 
of harmonic terms.* 

The laboratory at Columbia University is the first capable 
of performing intricate scientific calculations by machinery 
without any reading or writing of figures. The individual 
machines have highly specialized functions, but since each is 
capable of reading numbers from the punched cards and re- 
cording the results in the same manner it is possible to com- 
plete a long sequence of operations one after the other without 
writing down a single number by hand. 

A few of the specific astronomical problems which have 
been solved in this laboratory will be described in detail in 
later chapters, but to give some idea of their nature brief 
mention will be made of them here in general terms. 

One of these problems is that of planetary motion. It 
consists of the numerical integration of a set of three simul- 
taneous differential equations each of the second order. The 
solution is carried out with ten significant figures and at the 
rate of twelve epochs per hour. 

Tables of functions are constructed from their differences 
with great efficiency, either as printed tables or as a file of 
punched cards. From the latter type the values of the 
function may be taken automatically at the rate of many 
thousands per day. Calculations involving such tabulated 
functions are efficiently performed. 

The manipulation of harmonic series is useful in many 
kinds of theoretical work. Such operations consist of numer- 
ical harmonic analysis and synthesis and the multiplication 
of harmonic series. 

The formation and the solution of normal equations, the 
computation of correlation coefficients, of averages, etc., are 
natural applications of the punched card method. 

* On the Construction of Tables by Interpolation. L.J. Comrie. Monthly 


Notices, R.AS., Apr. 1928. The application of the Hollerith Tabulating Machine 


to Brown's Tables of the Moon. L. J. Comrie. Monthly Notices, R.A.S., 
May 1932, } 


INTRODUCTION 


In one program the punched card has been found to be the 


best means of recording astronomical observations; observa- 


tions recorded in this manner may be automatically reduced. 


The results of punched card calculations are easily listed 
The tabulator record may 


for reference or for publication. 
be used not only for printer’s copy but may in many cases be 


reproduced directly by photography. 


CHAPTER II 
THE PUNCHED CARD AND THE MACHINES 


THE OPERATING PRINCIPLE OF THE HoLLERITH MEtHOp 


5555559559999 99009 


In the Hollerith method of computation the numbers are 
recorded by means of holes punched in cards, and are read 
into the machines by means of electrical contacts made 
through the holes. In all the machines the principle is the 
same. The tabulating card, acting as an insulator, passes 


between a wire brush and a brass roller. 
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Fic. 1. Reading Brush of Punched Card Machine 


A hole punched in the card permits the brush and the roller 
to make contact, thus completing an electric circuit which 
contains an electromagnet. The position of the hole *n the 
card determines the time when the magnet is operated, and 
all parts of the machine which are synchronized with the 
passage of the card will therefore perf em different functions 
depending upon the position of the hole in the card. 
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In the sorting machine this magnet opens a chute along 
which the card slides until it reaches the proper receptacle 
In the adding machines, the digits on the card are added 
In the multiplying machine the multiplicand and multiplier 
are read from the card. 


?. THe TABULATING CARD 


The tabulating card (fig. 2) has 80 vertical columns with 
12 punching positions in each column. Of these, ten are 
indicated by the printed digits 0 to 9 which correspond to the 
digits of the numerical data to be punched. The 11th an 
12th punching positions at the top of the card (not indicated 
in fig. 2) are commonly referred to as the ae ead ty" 
positions. The “x” punch is usually used to control the 
operation of the various machines. For example, negative 
sere sth differentiated from positive ones by means of 
the “x”. ‘ 
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Fic. 3. Tabulating Card Partially Punched 


a Figure 3 shows a card with data punched in the first six 
line ste The punching reads:1 3.8 y 6 «90. The vertical 
lines have nothing to do with the punching, and are described 
in the next paragraph. 


. F . 
or the alphabetic machi 
Z machines ea TO AES SEL ‘ 
see sewn ne each letter is indicated by two punched holes 
’ r “y” being combined with one of the digits 0 to 9 


See sec. 4, 


Tue PuNCHED CARD AND THE MACHINES 7 


The first step in the use of a card for a particular record 
is the designation of groups of columns as “fields.” Each 
field defines a section of the card in which one particular type 
of information will always appear. If a given problem is 
large enough to warrant the expense, the cards may be ruled 
at the factory to indicate and label the fields.* Most of the 
computation of the Thomas J. Watson Astronomical Com- 
puting Bureau is done on a card with standard ruling designed 
only to assist the eye in reading from desired columns. 
Figure 3 shows the standard tabulating card used by the 
Bureau. The 80 columns are divided into fields of ten, with 
further subdivisions of three and one. 
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Fic. 4. Tabulating Card with Interpretation and Special Ruling 


Figure 4 shows a card which was ruled for an extensive 
project. There are over 30.000 cards in the set and ten 
duplicate sets are to be used at different places so the extra 
expense of ruling was justified. The numbers at the top are 
not part of the ruling but were added later by the ‘*inter- 
preter ’’ (see sec. 3). 


3. Tue MAacuINes OF THE Tuomas J. WATSON 
ASTRONOMICAL COMPUTING BUREAU 
The machines of the Bureau are, with one exception, 
standard accounting machines. Where necessary, circuit 
* There is no charge for the ruling but the preparation of a new electroplate 


for a change in ruling costs about $25. 
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changes were introduced, but the original functions of the 
machines have not been disturbed. The exception jis es, 
“Calculation Control Switch’’ which controls the hime 
of several other machines in the complicated process of ad 
ical integration. In describing these machines we shall ices 
fore give the trade name of each machine and note ve 
changes that have been made. ‘ 


The Card Punch 
The card punch is used to transfer the given data of a 
problem from the manuscript to the punched cards. The 3 
a . = ah - ONE 
used for most of the work of the Bureau is a ‘‘ motor-driven 
electric duplicating key punch.” In this model the cards ar 
ards are 


Duplicating 
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Keyboard 


Fic. 5. Duplicating Ke. 
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The operation of the punch is much like that of an ordi- 
nary typewriter. The card is fed into the machine with the 
first column under the punches. The depression of one of 
the twelve punching keys causes the corresponding hole to 
be punched in that column, and the card to move to the next. 
Columns may be skipped by pressing the spacer key or the 
card entirely ejected by pressing the eject key. Provision is 
also made for automatic spacing and ejecting on any pre- 
determined column. 

If a prepunched “ master card” is placed in the duplicat- 
ing rack, the information from the master card will be copied 
onto the card in the punch. Thus, for example, we might wish 
to punch common information in columns 1-5, 13, 14 of a 
large number of cards, columns 6-12 being different for each 
card. The common information is punched in the * master 
card” and it is placed in the duplicating rack. As soon as a 
card is fed into the punch, the duplicator copies columns 1-5, 
the operator punches columns 6-12. and the duplicator c ypies 
columns 13,14. When part of a card is spoiled the duplicator 
is used to copy the go yd portion. 

The speed of the punching operation depends not only 
upon the operator but upon the legibility of the manuscript, 
and the amount of coding required. The duplicator operates 
at the rate of ten columns per second, and about one and a half 
seconds are required to feedacard. A general idea of punch- 
ing speed may be obtained from a recent job where about 
half of each card was filled, part by the duplicator, and with 
a small part of the ‘nformation in code. On 30,000 cards'an 
average of 200 per hour was maintained. On another ex- 
tensive job, where about half the card was punched but con- 
siderable coding and rounding of numbers was required, the 
operators averaged less than 150. For straight numerical 
punching from go 1d typed manuscript 125 eighty-column 
cards should be punched per hour. 

The Verifier 

The mechanical verifier is used to check the punching of 

the initial data on the cards. It is operated in the same 
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manner as the key punch. 
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punching positions on the card and the thirteenth called the 
“reject” for cards which had not been punched in the column 
of the sorting. 

Provision is made for closing off any of these receptacles. 
For instance, if we wish to pick out only the cards with an “x” 
punch in the required column, we can shut off all the others 
so the cards without an “x will go into the reject. 


Fic. 7. The Sorting Brush in Operation 


The counting attachment gives the number of cards which 
goes into each compartment, the total for the group, and a 
grand total of all the groups. Switches are provided so that 
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counting may be done without sorting, sorting without count- 


ing, or sorting and counting. 
Figure 7 shows the movable sorting brush which may be 


quickly set on any desired column. | It also shows the card 
in the position where the actual sorting takes place. At this 
point the card passes into a chute leading to the proper 
receptacle. The arrow shows the brush making contact 


through the hole in the card. . 
Sorting is done at the rate of 20,000 cards per hour. 


The Tabulator 

The tabulator is essentially a high-speed adding machine 
of large capacity which reads the numbers to be added from 
holes in punched cards. When equipped with a printing 
mechanism it is called an accounting machine. ‘This printing 
mechanism will print information from each card or the accu- 
mulated totals only. 

The feeding mechanism of the tabulator has eighty brushes 
of the type shown in figure 1, one for each column of the card. 
As the card passes these ‘‘adding” brushes the adding wheels 
in the “counters” revolve and the electromagnets contre sling 
the addition are actuated by the electric contacts made 
through the holes in the cards. 

The electric circuits from the adding brushes to the 
electromagnets of the adding mechanism pass through a plug- 
board similar to that of a telephone exchange so that it is 
possible, by means of plugwires, to select the positions on the 
cards from which the numbers are to be taken and the posi- 
tions in the adding mechanism where they are to be added. 

The machine used by the Bureau is the so-called ‘type 
3S Direct Subtraction Accounting Machine.” It has been 
modified to some extent. The standard machine of this type 
has five adding units or “counters” each with eight-digit 
capacity, whereas this one has four ten-digit ones. A few 
special circuits have been added for the application described 
in chapter XII, but these do not interfere with the normal 
operation of the machine. It has the following features W hich 
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are optional on all standard machines: “automatic plug- 
board,” ‘major and minor control,” and two “class selectors.” 
The automatic plugboard is an arrangement whereby the 
plugboard may be removed from the machine for convenience 
in wiring. Several extra plugboards are provided so that a 
number of wiring arrangements may be quickly available. 
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Fic. 8. Subtraction Tabulator 


Automatic “‘control’’ is a feature of most tabulators. For 
control purposes the feeding mechanism contains, in addition 
to the “adding brushes” mentioned above, a second set of 
eighty brushes known as the ‘‘control brushes.” While one 
card passes the adding brushes the following card passes the 
control brushes. During the feeding process it is therefore 
possible to compare the desired part of each card column by 
column with that of the preceding. The principal use of the 
control is to enable the machine to sense when one group of 
cards has been finished and the next is about to begin. For 
example we might have cards numbered consecutively from 
one to one thousand and wish to tabulate them in groups o} ten. 
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We would wire the controlling mechanism for the columns 
on the card containing the tens and hundreds digits. In the 
two columns selected for the cont rol the first nine cards would 
match and would ‘pass control.” The tenth card however 
would differ in the tens position and would therefore ‘‘ break 
control” and start a new “ control group.”” This break in con- 
trol can be used to cause the tabulator to perform various 
functions such as printing totals from the counters and clear- 
ing them for the next group. By the use of major and minor 
control it is possible in the above example to tabulate by tens 
and also by hundreds. Any number of columns up to sixteen 
may be used for control. 
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Fic. 10. Feeding Mechanism f Tabulator 


Figure 10 shows the card feeding mechanism of the tabu- 
lator. The cards to be tabulated are placed in the hopper at 
the top. They are fed one at a time from the front of the 
pack and moved downward past the control and adding 
brushes, and are finally stacked in a second hopper. 
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The class selector makes it possible to add in a given 


} : Lif t fields on the card, dependi purpose is to punch on a new card the totals accumulated in 
ro different neids > Card, depending 
{ counter from one of two differ : 


the counters of the tabulator. It may be operated by a 


om a fartal . 
| upon whether or not the card _ ef wai i : ie irae major or a minor break in control. A plugboard on the 
i The use of two class selectors in series makes it possible to summary punch makes it possible to select the positions from 
I select from three fields. De ca he mad: the counters to be punched in any column of the card. By 1 | 
iF The individual counters oe ta es e a “ae Ate sls throwing a switch on the summary punch it is possible to | 
ii add or subtract the numbers read from a Cha Geeending Bike lM change the signs of all the numbers punched from a given ! 
i whether or not the card contains an “x” in a given column, h onpnter: 
| i The “x” controls are wired through the plugboard and hence i Information common to many cards may be punched from 
| | give great flexibility. For instance we mai add the “x cards i a master card in the duplicating rack as with the key punch, 
; Hl) and subtract the “no x’’ ones, or add the “x cards and omit so that part of the information on a given “summary card” 
im the “nox.” By the use of “y's” in two columns it is possible i may come through the tabulator and the rest from the master | 
fil to subtract algebraically the numbers on one set of cards from é card, 
i= those on another (see app. B). In addition to reading the hy The summary punch operates at the rate of ten columns i 
Ai numbers from the cards a special circuit is provided for reading per second, and 1.5 seconds are required to feed a new card. | 
HH a “1” with the passage of each card. When this circuit is ; : ani 
‘e wired to a counter it gives a count of the cards tabulated. if tke simaacaanal 
|| The printing mechanism consists of five printing “banks”: F [he multiplying punch reads two numbers from a card, 


one connected to each of the four counters and one, known asa F forms their product and punches it on the same card. The 
“list bank,” which prints directly from the cards without ! crossfooting’’* feature makes it possible to add or subtract 
i adding. It is thus possible to print fifty digits from each as well as to multiply two numbers from a card. The two 
card and to add forty of them in the counters simultaneously. machines used by the Bureau have the crossfooting feature, 
When each card is listed separately the process is known as and OMe of them has an automatic plugboard. 

“listing” but when only major or minor totals are printed lhe following operations are possible: 


the term “‘tabulating”’ is employed. A eB ASB CED ie 


In connection with the use of control groups there is a AM BLAC ASRS OD 
process known as “group indication.’” This means that for At BoC 


each control group the numbers selected for group indica- 
tion will be printed or added only once. Information from 
fields other than those used for the control may also be wired 
through the “group indicator.”’ 


The possibilities of the multiplier are best understood a 
from a description of its construction and operation. 
The feeding mechanism is similar to that of the tabulator. li, 


9.000 is per hour There are 80 reading brushes which correspond to the adding 
Pe eaten the rate of 7,000 cards per ho brushes of the tabulator and which are wired to the plug- 


ne ; 2 es eae waciitred of 
wie tantating, and 4,500 while listing. The time requ board. The eighty control brushes of the tabulator are here 
to break control and reset is 1.1 seconds. replaced by three or four movable brushes which may be set 
| The Summary Punch : * Crossfoot is an accounting term which means to add nimbers across a line 
The summary punch is a motor-driven duplicating key rather than down a column. 
punch attached by an electric cable to the tabulator. Its 
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in the desired columns. Their function is limited to the de- 
tection of “x” punches, and they are commonly referred to as 
“x” brushes. When a card has passed the brushes it goes 
into a card punch where the answer is recorded on it. 

The multiplying machine contains seven counters of the 
type used in the tabulator: six of eight-digit capacity and 
one of ten. They are referred to as (1) multiplier; (2) multi- 
plicand; (3)-(4) right-hand components; (5)-(6) left-hand 
components; (7) summary products. Counters (3)-(4) are 
connected in series to form one sixteen-digit counter, and the 
same is true of counters (5)-(6). Counter (7) is the ten-digit 
one. 

In straight multiplication the operation is as follows: As 
the card is fed past the adding brushes, the multiplier and 
multiplicand are read from it into counters (1) and (2) respec- 
tively. The multiplication then starts. The first digit of 
counter (1) multiplies each digit of counter (2) electrically 
and at thesame time. Thus if the first digit of the multiplier 
is a 9 and the number in the multiplicand is 7638, the 
partial products 9 X 7 = 63, 9 X 6 = 54, 9 X 3 = 27, and 
9 x 8 = 72 are formed from the “electrical multiplication 
table” all at one time. These partial products together with 
those which will come from the remaining digits of the multi- 
plier must be accumulated in the 16-digit counters of the right 
and left-hand components. The reason for two sets of coun- 
ters is that only one digit at a time can be added on a given 
adding wheel. The right-hand components of the above 
operation are of course 3, 4, 7, 2 and the left-hand components 
are 6, 5, 2, 7. 

After the multiplication by all the digits of the multiplier 
is completed, the contents of the right-hand components 
counters are transferred to the left-hand components counters. 
The results are then punched on the card and also transferred 
to the summary products counter. This counter accumulates 
the sum of all the products of the group. 

In punching a given product on the card, the entire 
product is usually not required. For instance in forming the 
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product of two eight-digit numbers we usually want the prod- 
uct to eight figures. In dropping the extra digits the last 
figure is usually “rounded” to the nearest unit. This round- 
ng is automatically done by the machine with the use of the 
‘“‘half pick up”’ whereby a ‘‘five’”’ is automatically added in 
any desired column of counter (4) immediately after the card 
passes the reading brushes. 

The crossfooting operations are accomplished by reading 
into counters (4) and (6) at the same time as (1) and (2). 
When the numbers are being transferred from (5)-—(6) to 
(3)—(4) they may be transferred as true figures or as comple- 
ments, thus giving the subtraction crossfooting. On some 
of the later models it is possible to read also into counter (7). 

For purposes of checking, a switch is provided which 
interchanges the multiplier and multiplicand, and spaces the 
card without punching. The sum of the products is then 
found for comparison with the value obtained during the 
previous run. 

In ordinary multiplication or crossfooting all the counters 
except the summary products counter are cleared after each 
card. When a common multiplier is to be used for several 
cards it is possible to retain the multiplier in the counter. 
This may be done in two ways: “fixed multiplier,” or “‘rate 
card.”” In the case of the fixed multiplier the multiplier 
counter must be reset by the operator. In the “rate card” 
method each card containing a new multiplier is identified by 
an “‘x’’ punch. This “‘x’’ punch causes the old number to 
be cleared out and a new one to be put into its place. The 
multiplicand from each card is thus operated upon by the 
multiplier from the last preceding rate card. 

Special features of our machine, aside from those men- 
tioned in chapter XII, are: 

(1) A device for reading the algebraic signs of multiplier 
and multiplicand from ‘“x’’ punches on the card, and indicating 
the sign of the product in the same manner. 

(2) Aspecial circuit for forming the quantity A + (B X C) 


: 
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with the sign depending upon whether or not the card carries 
an “x"’ punch in a given column. 

(3) Provision for additional crossfooting. The standard 
machine does not provide for the use of all sixteen possible 
positions in counters (3)-(4) and in (5)-(6). 

The machine operates at 150 cycles per minute. Two 
cycles are required to feed the card, one for each digit of the 
multiplier (zeros excepted), and one for the transfer from 
counters (3)-(4) to (5)-(6). The speed is independent of the 
number of digits in the multiplicand or in the crossfooting. 
Punching requires one second for ten digits, but part of the 
punching is done on the feed cycle of the next card. Thus 
eight by eight digit multiplication is done at the rate of 
730 cards per hour when ten digits of the answer are punched, 
and at the rate of 870 for “non-punch.”” Three by eight 
multiplication, “‘non-punch,” is performed at the rate of 1,350. 


The High-Speed Reproducer 

This machine is used to transfer information from one 
card to another. It is much faster and more flexible than 
the duplicator of the key punch. An automatic plugboard 
allows convenient changes of wiring. Comparing relays check 
the new card against the original during the same run of the 
cards. 

A B 


READING COMPARING 
BRUSHES BRUSHES 
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Fic. 11. Schematic Diagram of High-Speed Reproducer 


There are two complete card feeding assemblies similar 
to that of the tabulator. In one of these feeds the control 
brushes have been replaced by eighty punching knives. Fig- 
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ure 11 shows the paths of the cards through the two feeding 
mechanisms. The two may be operated independently, or in 
synchronism with a card in feed “B” in exactly the same 
phase as the corresponding card in feed ‘‘A.” In addition 
to the brushes shown in the above diagram there are several 
movable ‘‘x” brushes of the type used in the multiplier. 

A file of cards may be reproduced by passing it through 
feed “A” while new cards are passing through feed “B.” 
The reading brushes of “‘A”’ are wired to the magnets operat- 
ing the punches of ‘“‘B,” and the comparing brushes of “A” 
and “B” are wired to the “comparing relays.” The purpose 
of the comparing relays is to compare the holes in the two 
cards which are passing the brushes connected with the relay. 
If the holes do not match, the machine stops. 

“Gang punching” gives a large number of duplicates 
from one master card. Gang punching is done in feed ieee 
The master card is put in front of the new cards and the 
information from each card is passed back to the one following. 
For this operation the comparing brushes of feed “B”’ are 
wired to the punches. 

While feed ““B” is being used for gang punching, feed 
‘‘4”’ may be used independently for checking. 

Where several master cards are interspersed among the 
“detail” cards, an “x” punch in the master card will prevent 
its being punched. 

The reproducer operates at the rate of 6,000 cards per hour. 


The Interpreter 

The interpreter prints in type at the top of the card the 
information punched in the card (see fig. 4). A plugboard 
permits the selection of the columns to be interpreted and the 
position of the printed characters on the card. Interpreting 
may be done on the front or the back of the card. Our 
machine will interpret as many as 45 digits. The newer 
models will interpret more columns, but have smaller type. 

The capacity of the machine is 6,000 cards per hour. 
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The Calculation Control Switch 

The machines described above are all designed for com- 
putation where each operation is done on many cards before 
the next operation is begun. For the process of numerical 
integration the different arithmetical operations must be done 
in succession for the first step before proceeding to the next, 
For this work it must be possible to change quickly from one 
complicated setup to another. 

The calculation control switch controls the operation of 
the multiplier, tabulator, and summary punch so that a cycle 
of different arithmetical operations can be performed in rapid 


succession. 


4 Orner Macuines AND Devices Nor Now IN Use 
AT THE BUREAU 


Alphabetic Equipment 
For commercial accounting there is a complete line of 
punches, tabulators, and interpreters, which handle literal 
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Fic. 12. Alphabetic Tabulator Record 


as well as numerical data. Figure 12 shows information 


hc from punched cards by means of an alphabetic tabu- 
tor. 
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This method of recording literal information on cards 
and printing partial lists from the cards is an efficient means 
of furnishing current information from a file which is changing. 
For instance in the U. S. Coast and Geodetic Survey the 
descriptions of bench marks and their elevations could be 
kept in such a file. Of course where the information is pre- 
dominantly numerical it is possible to code a small amount of 
literal data and use the numerical machines. 

The alphabetic tabulator has many advantages over the 
numerical tabulator for numerical work. For instance some 
of the ‘‘carry-over” mechanisms of the counters are wired 
through the plugboard so that the counter capacity may be 
arranged to suit the problem. The chief disadvantage is the 
increased cost. 


The Collator 


The collator is a new type of machine which is important 
for scientific work. It is really a mechanical filing clerk 
which may be used to insert or remove cards from a file. 
For example you have a file of 10,000 cards numbered con- 
secutively in columns 1-5. You have another set of “detail” 
cards containing any 2,000 of the numbers from 1-10,000. 
The collator will match each one of the detail cards with the 
proper file card and remove it from the file. The collator 
will later insert each of the file cards in its proper place. 

In the work of the Bureau the sorter has been used to 
perform many of the functions of the collator, but of course 
with less efficiency. 

The Test Scoring Machine 

This machine is not a punched card machine, and it is not 
a computing machine, but it is an interesting device which 
might well have uses in science. 

It was designed to score the so-called “new type” exami- 
nations. Each question in the examination is answered by 
placing a lead pencil mark in one of five positions. The com- 
pleted test sheet will thus have a pencil mark for each of the 
100 questions. The test scoring machine compares these 
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ks with similar ones on the master sheet and enstan- 
mar usly gives the total score on a voltmeter. rhe score 
cee be computed in several ways with various weights as- 
‘oned to the different sections of the examination. | 
sign ration of the machine depends upon the electric 
ie ay a lead pencil mark rather than upon contact 


made through a hole in a card. 


Miscellaneous Devices . 

With a high-speed listing device a standard tabulator will 
list 6,000 cards per hour instead of 4,500. Bile. 

A tabulator may be equipped with special fractional 
wheels” for adding fractions or degrees, minutes, and seconds. 

We have seen above how the operation of the counters in 
the tabulator may be controlled by an “x” punch in a given 
column of the card. The “digit emitter”’ is a device which 
will operate the counters according to the digits punched in 
a given column. For example if our tabulator had a digit 
emitter we could add in counter (1) information from all the 
cards with a “2” in column 17, in counter (2) from those with 
a “7” in column 17, etc. This device lends itself to a great 
many useful tricks. 

Intermediate control provides an additional subclass in 
the tabulation. 

The collating device on a sorter will deal cards for a bridge 
game or similar purpose. 

The high-speed summary punch consists of a gang punch 
attached to the tabulator instead of a key punch, It is 
much faster, but of course more expensive. 


PART II]. MACHINE METHODS 
CHAPTER III 


GENERAL CONSIDERATIONS OF THE 
PUNCHED CARD TECHNIQUE * 


1. INTRODUCTION 


The successful and efficient application of the machines 
to a scientific problem requires careful consideration of human 
and mechanical elements. The machines are not infallible, 
so numerical checks must be provided. The machines are 
not infinitely efficient, so it is necessary to take advantage of 
all short cuts. The obvious direct method is frequently much 
more costly than necessary. It is often possible to substitute 
a little manipulation of the cards by the operator for long 
machine processes. On the other hand the operator should 
not be assigned tasks which are psychologically difficult to 
perform without error. Part III of this book illustrates these 
features, but a few of them may be summarized here in general 
terms. The planning of an extensive program thus requires 
a careful analysis of many factors in the light of all available 
data and experience. 

It should be remembered that electric computing machines 
of the keyboard variety are now so efficient that unless each 
number required in the computation is to be used several 
times the punched card method would not be efficient. The 
advantage of the punched card method increases with the 
number of uses to be made of each number concerned. 


2. NUMERICAL CHECKS 
The calculation must be so arranged that a check is pro- 
vided at each step. The machines are remarkable in their 
* See also Practical Applications of the Punched Card Method in Colleges and 
Universities, edited by G. W. Baehne. Columbia University Press, 1935. 
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accuracy, but errors do occur. — It must be remembered that 
any failure in the entire machine may in a short time ruin 
large quantities of work. It is not sufficient to have a check 
at the end of a problem requiring weeks to complete, but each 
operation must be checked before proceeding to the next, 

While it is not the purpose of this book to discuss tabulat- 
nique which is part of the stock in trade of the average 
operator, it is not out of place to mention some things which 
are more important in scientific than in most other work with 
the punched cards. In scientific computation a given card 
may receive more varied treatment than in other work. 
More different machines are used, the card may be run 
through the machine with the back side forward, etc. For 
these reasons the cards must be treated with great care and 
the adjustment of the machines must be carefully checked. 
Each new set of cards should be checked on the “card gauge” 
to see that they are the exact size. The cards should be 
carefully stored in files which hold them properly. Above 
all the operator must insist that the machines be kept in 
perfect adjustment so as to prevent “off punching.” This 
can be avoided if cards coming from all the machines are 
frequently tested with the gauge. These are all standard 
techniques but carelessness is fatal to successful work of the 
kind described here. 

The first thing to be checked is naturally the punching of 
the original cards from the manuscript. Three methods are 
available for this: the use of the mechanical verifier, proof- 
reading, and the use of check sums; which is to be used 
depends somewhat on the nature of the problem. [If the 
initial data are not extensive and the problem is long and 
costly all three methods may be employed. If the data are 
extensive and of a rough nature such as are used in many 
statistical discussions there is no point in trying to find errors 
which would not materially affect the result. 

The mechanical verifier in the hands of a conscientious 
operator is probably the best means of checking ordinary 
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numerical data. It does not, however, detect double punch- 
ing in a given column. 

The cards may be listed on the tabulator and proofread 
against the manuscript. This has two advantages. The 
comparison may be made by someone who is familiar with 
the data but is not a skilled key punch operator. Again, if 
the original record is not written out in full, or if the informa- 
tion about each case is not complete, the method of proof- 
reading is simpler. 

Checks which depend upon sums or upon sorting are 
desirable because of the efficiency with which these operations 
may be performed on the cards. Check sums detect mic- 
placed cards and errors of transposition. Another advantage 
of check sums is that they may be carried on to the next step 
in the work. 

Ordinary sorting is checked by looking through the holes 
for small groups, and by passing a sorting needle through the 
larger ones. This checks each sorting as it is made, but does 
not insure against getting them out of order later. For 
example in putting the cards in order on several digits, sorting 
is done first on units, then on tens, etc. While sorting on the 
tens it is essential that the units keep their proper place. If 
the stacking mechanism of the sorter is out of adjustment one 
card may stand on edge long enough for the next one to pass 
under it. 

There is no direct infallible check on the listed numbers or 
printed totals from the tabulator, but with proper precau- 
tions the probability of such errors can be made very small 
indeed. While listing with the tabulator, the numbers being 
listed should also be added in the counters. The totals accu- 
mulated in the counters may then be compared with similar 
totals obtained by the tabulator. This insures that the num- 
bers were properly read from the cards. The failure of the 
printing mechanism itself is rare and usually of such a nature 
that it is easily detected. 

In tabulation it is desirable to punch the totals with the 
summary punch. The addition of these summary cards then 
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perfect check on the tabulation. It is sometimes 


wee k done on the tabulator to vary the 


‘rable in checking wor Rees 
paar so as to eliminate the possibility of the tabulator 
Pp 


king the same error twice. For instance listing and tabu- 
mak 


lating are done at different speeds of the machine. The 
ating erchanged. The counters may be wired 


int 
ig oy at of all the numbers. The cards may be 
reversed from left to right so as to use different adding brushes. 
The number of such precautions to be taken depends upon 
how serious an undetected error would be. 

In checking summary cards, most totals may be checked 
by straight tabulation. The parts, however, which come from 
the group indicator require a special technique. The plug- 
board is so wired that the tabulator does not break control 
on each control group, but the counter which is adding the 
group indication operates only from the impulses which in 
the summary punching caused the break in the control. The 
wiring for this operation is given in appendix B. 

In transferring information with the high speed reproducer, 
the comparing brushes should be used to check the cards as 
they are made. Tabulation of the cards is, however, a desir- 
able additional safeguard. It is possible for a card to feed 
incorrectly into the reproducer, to be incorrectly punched, and 
still to pass the comparing brushes. 

Ordinary multiplication is checked by running the cards 
through again with the machine set for ‘‘check.”” As stated 
in chapter II, section 3, this interchanges the numbers in the 
multiplier and multiplicand counters, and accumulates the 
sum of the products. This form of check is practically perfect 
as far as the multiplication is concerned. It leaves the slight 
possibility of an error in the transfer from right-hand com- 
ponents to left. In the newer multiplying machines the 
product formed in the check multiplication is compared with 
the punched value on the card. This checks the individual 
cards and locates the error at once, whereas in the older 
method the fact that the sum of the products does not agree 
simply tells that an error has been made. 
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When many products are formed with a constant multi- 
plier they should not be checked as above. The sum of the 
products should equal the product of the constant multiplier 
by the sum of the multiplicands. Usually not all the digits 
of the product are punched on the card so the check on the 
tabulator is not exact. However, when the multiplying is 
being done additional digits of each product may be accumu- 
lated in the summary products counter to check the last 
required digit. Suppose for example you have a group of 
500 cards with the multiplier and the multiplicand each of 
eight digits and the products punched to ten. The summary 
products counter has ten-digit capacity so we must drop the 
three largest digits in order to leave room for the three addi- 
tional smaller ones. The larger digits in the sum of the 
products are of course obtained by tabulating the punched 
products. 

Where many small products not involving a constant 
multiplier are to be checked we can form the sum of the 
products by the method of section 5. When this method is 
to be employed the cards should be sorted in order of the size 
of the multipliers before multiplication. Small groups may 
then be checked on the tabulator as they come from the 
multiplier. 


3. SuHort Cuts IN MULTIPLICATION 
Multiplication is one of the most expensive operations, as 
may be seen from the operation times given for each machine. 
It is therefore desirable to find means of avoiding it whenever 
possible. When the individual products are not required but 
only their sum, the method of section 5 should be considered. 
When the multiplying punch is used it should be used to 
capacity. For instance, it is possible to do more than one 
multiplication at a time if the numbers are small. It is only 
necessary to avoid one product overlapping the other. Asan 
example consider the formation from a, } of a®, 2ab, b*, where 
a, b are two digits each. Suppose for illustration @ = 96, 
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b = 84. The multiplication is as follows 
a b 


a > 
96000084 Multiplicand 
96000084 Multiplier 


i al 
384000336 


768000672 
576000504 
864000756 i 
9216016128007056 Product 
a 


Unfortunately this method will not work for two three-digit 
numbers. It will work for a two and a three. Where a 
great many such products are to be formed and a and 6 are 
numbers each with not more than three digits it is possible 
to sort out before multiplication those for which a is less than 
three, then those for which a is 3 but b is less than three. 
The remainder for which both a and 8 are three digits must 
be done the long way. This technique may also be used as 
follows: a + 6 in multiplicand and c + d in multiplier giving 
ac, ad + bc, bd in the product. 

It should be remembered that when a product is being 
formed the addition of a crossfooting operation usually 
requires no additional machine cycles. The two operations 
should therefore be combined whenever possible. It is pos- 
‘sible to perform for example a + 5, ¢ X d at the same time. 

Where large numbers of small products are required they 
may be obtained with the sorter and gang punch. The cards 
are sorted into groups according to the multiplicand and 
multiplier. A small table gives the product for each group 
and it is put on the cards with the gang punch. This method 
pays only when most of the groups contain a considerable 
number of cards. 
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4, SpectIAL SorTING TECHNIQUE 


For many calculations it is necessary to sort the cards 
into groups with characteristics depending indirectly upon 
the information punched in a given field. As an example of 
such sorting consider the problem of sorting m cards, numbered 
consecutively, into » groups according to the remainder left 
after dividing the number by m. The general method requires 
as many sortings as there are digits in m. Let us illustrate 
with the case where there are 9,999 consecutive numbers and 
the period m is 13. The four digits of the numbers are 
punched in columns 1-4. Sort first on column one, forming ten 
packs of cards. Then sort each pack on column two. The 
first pack when sorted on column two gives 00, 01, 02, --- 09. 
The second gives 10, 11, --- 19. When we come to pack 13 
we combine it with pack 00, combine 14 with 01, and so on. 
The sorting is now done on column three in the same manner, 
and finally column four. The final result gives 13 packs of 
cards. The first contains cards 00, 13, 26, ---, 9997, the 
second 01, 14, 27, ---, 9998, and so on. 

The general method is not necessary for all periods. For 
periods 2, 5, 10, only one sorting is required; for periods 4, 20, 
25, only two; and for 8, only three. 


5. THe FoRMATION OF SUMS OF PRODUCTS WITHOUT 
THE INDIVIDUAL PRODUCTS 


In the method of least squares, and in the computation 
of correlations we frequently require the sum of a great many 
products, the individual products being unnecessary. Such 
sums of products may be very efficiently formed with the aid 
of the sorter and tabulator.* 

The process is known as ‘progressive digiting”’ and may 
best be explained by a simple example. 


*The Mendenhall-Warren-Hollerith Correlation Method. Columbia Uni- 
versity Statistical Bureau Document No. 1, Also The Use of Machine Factoring 
in Multiple Correlation, A, E. Brandt, Journal of the American Statistical Asso- 
ciation, September, 1928. 
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We require the sum of the following products 


ix2= 2 
ixs= 8 
2x3= 6 
ax1l= 2 
2x5= 10 
3x9= 27 
5x6= 30 
5x2= 10 
6x7= 42 


ol 
Sum = 137 

This sum may be found by addition only. Simply start at 

the bottom and add the seven six times, the two and the six 

five times, etc. The required additions are given below. 

The last column contains in each line the sum of the numbers 


in the adjacent column. In other words it is the ‘‘ progressive 
total” of the numbers in the preceding column. 
7= 7 
2+6= 15 
0= 15 
9=9 24 
54+14+3=9 


—_ 
Ww — WwW 
a | ww wo 


The sum of the numbers in the last column gives the required 
sum of the products by simple addition. 

The additions of the type outlined above are done with 
great efficiency on the tabulator. The cards are put in order 
with the sorter, and the automatic control causes the printing 
of a total after each group. Blank cards must of course be 
inserted for the missing groups. If the summary punch is 
available the progressive totals are punched on cards and 
added with the tabulator; if not they must be added from the 
printed tabulator record. 
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The great efficiency with the tabulator comes from the fact 
that the product of one number by several others may be 
treated simultaneously. For example, in the method of least 
squares or for correlations we have several columns of num- 
bers, say a, b,c, d,--- and we need the sum of products 
[aa], [ab], [ac], [ad]. These quantities are formed simul- 
taneously on the tabulator, [aa] in counter (1), [ab] in 
counter (2), [ac] in counter (3), and [ad] in counter (4). 
When the numbers in a, b, c, d are small, one counter could be 
used for a and 6, thus doubling the capacity of the machine. 

In the simple example above we have assumed a, b each 
of one digit. Two digit numbers could be treated in the 
same way with, at most, 99 groups; or the work could be done 
in two steps with 9 groups each. For this latter method the 
procedure is as above for the first step except that both digits 
of the multiplicand are added. When the first step is com- 
pleted the operation is repeated for the tens digit of the 
multiplier. The results of the two steps are of course com- 
bined by addition. This method may obviously be extended 
to three-digit numbers. 

The following example illustrates the two-digit case. 


Se S 

2x13 = 26 

52 f= 2) 

ll x 16 = 176 

12 3 ‘= 36 

23x 2 46 

Sum = 313 

For the unit digits we have 

2+ 7=9 9 
3 4+.13:= 16 25 
146+ 8=24 49 
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For the ten digits 
Omer. 2 2 
3+16=19 21 
23 


Combining, 10 X a4 4. 83 = 313. 


6. SpecIAL CARDS FOR BREAKING CONTROL 


In many problems it is desirable to form control groups 
according to information not directly punched on the cards. 
For instance we have on the detail cards a set of consecutive 
numbers and wish to form control groups of threeeach. This 
may be accomplished by sorting into the detail cards a special 
set of ‘control cards.”’ These control cards must be punched 
so they may be sorted into the detail cards and come at the 
desired places. In addition to the punching required for sort- 
ing the control cards into their proper places they must contain 
an “x” in some given column. This “‘x’’ when used in con- 
nection with a special wiring of the plugboard (see appendix B 
will cause only one break in control for each control card 
instead of the normal two. An example of the use of control 
cards will be given in chapter V, section 3. 


7. SUGGESTIONS FOR CLERICAL Work 


Where many computations are being made simultaneously, 
some of which extend over long periods of time, general 
schemes of identification and indexing are desirable. 

One method of identification is by means of colored cards. 
In our laboratory manila cards are used for most of the com- 
putations because they are the cheapest. Colored cards are 
used in special cases. For instance w hen a card jams and 
must be reproduced the reproduction is made on a salmon 
colored card. The color aids in checking to see that the torn 
card was properly reproduced and also shows up in the file 
as a card which is more likely to contain errors than the ¢ thers. 

For a very extensive project it is desirable to have spe ial 


ruling and printing on the cards. Small files may be identified 
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by ruling stripes on the top edges of the cards. The positions 
of these stripes can indicate the number of the file. With 
three stripes about five hundred files may be quickly identified. 
Phe numbers should of course refer to a description of the 
file. These stripes are easily applied to a trayful of cards 
with aruling pen. The marks are visible on the cards without 
opening the tray, and also on the individual cards. 


a, 
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a SE OF SPECIAI 
= CONSTRUCTION AND USE OF SPECIAL 
THE Fs OF TABULAR FUNCTIONS 


1. INTRODUCTION 


The tabulator furnishes such an efficient means of con- 
« © J 


tructing tables of tabular functions that every computational 
5 s «< 


problem must be examined to see to what extent special tables 
may be used in its solution. In some cases the construction 
of several manuscript tables to be used in the usual way may 
furnish the most economical solution of the problem. In 
others the special tables should be prepared in the form of a 
punched card catalogue to be used in connection with other 


machine work. z 
The interval adopted for the tabulation depends upon the 


number of times the table is to be used. In some extreme 
cases it may pay to tabulate so that interpolation is unneces- 
sary. In this case “critical” tables are convenient. The 
most common tables are those in which the interval of tabula- 
tion is chosen just small enough that linear interpolation with 
first differences is legitimate. Where comparatively few 
values of the function are to be used, the interval may be 
chosen so large that second and higher differences are required 
in the interpolation. In view of the efficiency with which 
the tables may be built up on the tabulator, smaller intervals 
will in general be taken than if the tables had to be con- 
structed by more laborious means. 


1. SUMMATION WITH THE TABULATOR 


The tabulator will automatically build a table from the 
differences, and print argument, function, and difference with 
one run of the machine. The method is extremely efficient 
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and is self-checking. Consider as an example the construction 


of a table of logarithms beginning at 5000. 


Arg Log Diff 
5000 69897 
g 
5001 69906 


5002 69914 
5003 69923 
5004 69932 
5005 69940 
5006 69949 


In building this table on the tabulator three counters are 
used; the first for the argument, the second for the function 
and the third for the difference. The first counter builds up 
the argument as a progressive total by adding “1” from each 
card, and the second builds the function by adding successive 
values of the difference. The third counter is used only for 
the printing of the individual differences and is cleared or 
reset after each line is printed. 

The differences 9, 8, --- are punched on a series of cards 
by means of the key punch (unless they are the result of a 
previous calculation with the punched card machine). With 
the gang punch a “1” is punched in each of these cards. 
We also prepare one card with the first value of the argument 
and of the function, in this case 5000, 69897. The tabulator 
is wired so that 5000 is added in the first counter and 69897 
in the second. Further the “1” from the difference cards 
must be wired to the unit position of counter (1), and the 
field containing the difference must be wired to counter (2) 
and also to counter (3). The automatic control is set to 
break control on each card. Counters (1) and (2) are set for 
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“total, non-reset,” and counter (3) for ‘‘minor total and 
' The tabulator record would appear as follows. 


reset.” 
5000 69897 
5001 69906 9 
5002 69914 8 
5003 69923 9 


The value of the difference which is printed on any given line 
is the one used in obtaining the value of the function on that 
ine from the preceding. It is possible to set the machine to 
print the difference on the odd line of the table if desired. 

In order to check such a table it is only necessary to 
check the first and the last values of the function and argu- 
ment and the values of the differences in column three. Thus 
each line contains ten digits, only one of which requires 
personal attention. With the subtraction tabulator we can 
take care of positive and negative values of the function and 
difference. 

If the table is to be used in the form of a punched card 
file we use the same process but punch summary cards to 
match the above printed record. There is one important 
difference, however, and that is that we want on a given card 
the argument, the function, and the value of the difference 
which is to be used in interpolating between the value of the 
function on the card and the value corresponding to the next 
value of the argument. This may be accomplished by build- 
ing the table from the bottom rather than from the top. The 
printed record is then inverted but the cards can be sorted 
into the normal order. 

In the above example the differences were taken from a 
table which had already been printed. In the construction 
of a new table the differences are obtained without first 
computing the values of the function. Methods of obtain- 
ing the differences required to build up a table with a smaller 
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tabular interval than that of a given table will be discussed 
in chapter V. 

The process described above may be applied successively 
to build a table from second or higher differences. In this 
case if the summary punch is not available it is necessary to 
punch the results of the first summation on cards in order to 
use them as the differences for the second process. Where the 
table is built from second differences counter (4) would be 
used to print the second difference in the final table. 

In printing the higher differences of a table the arrange- 
ment depends upon the interpolation formula to be used in 
taking values of the function from the table. For instance it 
may be desirable to use one of the well-known formulae 
requiring horizontal instead of diagonal differences. Before 
building the table of the function from the first difference it 
is possible to move all the values in the second difference 
column up or down the table by means of the high-speed 
reproducer. For example, suppose we have a file of cards 
each containing a first and a second difference and it is 
desired to transfer to each card the value of the second differ- 
ence from the preceding card. For this operation feed ““B”’ 
of the high speed reproducer is used. The plugwires from 
the brushes are so arranged that numbers read by the com- 
paring brushes are punched in a new field on the card. 
Checking may be done on the tabulator by adding the num- 
bers from the two fields. If the cards are run through the 
machine in reverse order we get on each card the following i 
instead of the preceding second difference. If the interpola- 
tion formula requires the sum of two adjacent second differ- 
ences, this sum may be obtained by crossfooting with the 
multiplying punch. 


3. THE Use or TABLES IN THE ForRM OF PUNCHED CARDS 

For many problems the ease of taking values of a function 
from a table by the punched card method is a very important 
factor. This is especially true where the previous and follow- 
ing steps in the calculation are done by means of punched 
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cards. The method is most efficient when a great many 
values of the function are han iain A oe a Hrs completely used before reproducing. 

Let us consider a Ee ctcent E escetion, ae : WV The collator is much more efficient than the sorter for 
each card cpap Ses etidetail xa ica . combining the file and detail cards. In addition, the question 
ence. There is also a large t ; t Saag 7 wn of matching the fields of the argument does not arise, for the 
taining an argument for which the ae is © ssi and plugboard permits sufficient flexibility in the wiring. 
an “x” punch in some given column. Suppose further that One run on the sorter is required to separate the file from 


demands that an attempt be made to have the previous card 


the function is tabulated in the table with three digits in the the detail cards. The column with the ‘‘x”’ is suitable for this. 
— and alain es ie cceunent on chs The interpolation of the function for the last three digits 
ve six. ne e { of the argument is now done on the detail cards. This i 
detail card which correspond to the three digits of the argu- ia ane on the multiplier as an Ae tisiaiee te — 
ment of the table cards must be in corresponding columns. Ax+zBXC. The special ‘“‘x” brush on our machine makes 
The first step is to sort the table cards and the detail it possible to do all interpolation of this form without separat- 
cards together in order of the three digits of the argument, ing the plus ones from the minus. 
and with the table card preceding the detail cards in each | Where interpolation with second differences is required 
case. Thus for instance if the argument of the table runs the procedure is as above except that two runs on the multi- 
from 1 to 999 we have, after sorting, table card 001 followed plier are in general required. The particular values of the 
by all the detail cards 001; then table car d 002 followed by first and second differences which are to go on the card depend 
detail cards 002, and so on. upon the particular interpolation formula to be used (see ch. 
The entire set of cards is now passed through feed ‘“B” of V). Stirling’s formula is most convenient for this purpose 
the high-speed reproducer. The wiring is arranged so that but requires smaller third differences for a given accuracy 
information from each card is transferred to the next, pro- than Bessel’s. Stirling’s formula for this purpose would be 
vided the card being punched has an “x” in the prescribed written: 


column. Thus the values of the function and of the difference ’ , ” 
‘ SE as = (8’_172 + 8/12) ) 
are punched on each detail card following it, but nothing is Pah se FSP 
punched on the table card. The machine may be wired to 
punch only the cards without an ‘“‘x,”’ in which case the “x” This requires two runs through the multiplier of the form 


should be on the table cards. The “x’’ on the detail cards is A = B X C provided each card contains n; fo, 60/2, and 
preferable because in case the machine fails to pick up an (8'12 + 5'12)/2. Bessel’s formula 


“x” the table card will not be spoiled. The comparing m= 1) 
relays are used to check the punching. The tabulator may fo =fot nin tn~ qe (80 + 41”) 
of course be used for additional checking. 

In planning a given computation provision should be made is not as simple to apply but might be used in special cases, 
that the field on the detail cards containing the digits of the for instance when the third differences are just too large for 
argument which match those of the file should be in the . Stirling’s formula. 
proper columns, either when the card is direct or reversed. In view of the methods here suggested for interpolation it 


If the cards must be reproduced to obtain this result, efficiency . will be seen that the considerations leading to a choice of the 
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tabular interval are different from those used in connection 
with hand work. For example in interpolating with first 
differences it is only necessary to secure an interval where it 
is legitimate to do so, the size of difference being unimportant 
to the machine. For example, Peters table of natural func- 
tions designed for ordinary use is given for each thousandth 
of a degree, whereas for the card catalogue we need only 
hundredths. 

In some cases the table should be tabulated for such a 
narrow interval that part of the interpolation may be done 
with the sorter and the gang punch. (See chapter IX.) 

A list of the tables in the form of card files available at 
the Bureau is given in appendix A. 


CHAPTER V 


INTERPOLATION, MECHANICAL QUADRATURE, 
AND ALLIED SUBJECTS 


1. GENERAL REMARKS 


It is not the purpose of this chapter to discuss all the 
available formulae for interpolation, numerical integration, 
and numerical differentiation, but rather to indicate how the 
punched card method may be applied to advantage. The 
reader is referred to the many books on the subject, to the 
collections of formulae, and to the tabulated values of the 
various coefficients required in the different formulae. The 
following short list is suggested as a beginning: 


(1) Steffensen, J. F. Interpolation. Williams and Wilkins. 
1927. 

Comrie, L. J. Interpolation and Allied Tables. Re- 
printed from the Nautical Almanac for 1937. His 
Majesty’s Stationery Office. 

Comrie, L. J. On the Construction of Tables by Inter- 
polation. Monthly Notices of the Royal Astronomical 
Society, volume 88, April 1928. 

Pearson, K. Tracts for Computers, No. II. On the 
Construction of Tables and on Interpolation, Part I, 
Uni-Variate Tables. London, Cambridge University 
Press, 1920. 

Bower, E. C. Systematic Subdivisions of Tables. Lick 
Observatory Bulletin 467. 1935. 


All the well-known methods of interpolation and me- 
chanical quadrature are equivalent to finding a polynomial of 
degree » which represents » + 1 consecutive values of the 
function, and then computing the required quantity from this 
polynomial. In the formulae of the Lagrangian type the result 
is expressed as a function of the ™ + 1 successive tabular 
values of the function; in the difference formulae one value 
of the function and the values of » successive differences are 
employed. 
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The justification of this process may be established for 
known functions which have continuous derivatives of a 
certain order.* For other functions the computer must be 
satisfied by the apparent convergence of the calculations 
themselves or from a knowledge of the general behavior of 
the functions considered. 

The choice between the two methods depends not only 
upon the problem itself but upon the computation which 
precedes and follows it. For example we might use a differ- 
ence method instead of the Lagrangian where the differences 
had been formed in the previous work. 

The principal disadvantage of the difference method is 
that the work of forming the differences may in some cases 
be greater than the complete interpolation by the Lagrangian 
method. On the other hand the differences form a check of 
previous work which should not be overlooked. Another 
advantage of the difference method is that the differences are 
smaller than the functions and hence the necessary multiplica- 
tion is easier. 

The relative merits of the various methods have been 
frequently discussed, for example see references (3) and (4) 
at the beginning of this chapter. In addition to the classical 
arguments should be included those which are introduced by 
the punched card method. By way of illustration, consider 
two examples from chapter XII where the decision was easy to 
make. 

For the numerical integration of the differential equations 
of planetary motion the difference method was chosen because 
the check afforded by the differences was urgently needed, 
and because two extra places in the computation were pos- 
sible. In this case a single error would invalidate weeks of 
subsequent computation. The differences checked not only 
the integration itself, but the perturbations which had been 
previously computed. In the second place, the multiplying 
punch has only eight digit capacity which would limit the 
Lagrangian methods to eight figures. The difference method 
gives the possibility of ten. 

* See page 2 of reference (1). 
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In the interpolation of the perturbations for the integra- 
tion, however, the Lagrangian method is used. This inter- 
polation is usually to halves and quarters instead of to tenths, 
so that the efficiency of the summation process is not so great 
as for tenths. The previous calculations did not require dif- 
ferences so it would have been necessary to compute them 
just for the interpolation. Finally the integration which is 
the next step in the calculation furnishes an automatic check 
not only of itself but of the perturbations as well. 


2. NOTATION 

In this and other chapters we shall use a uniform notation 
for tabular functions and differences. A given table consists 
of successive values of the argument and the corresponding 
tabular values of the function. The quantities formed by 
subtracting each successive value of the function from the 
following one are’ called the differences of the first order. 
The differences of the second order are formed by similar 
operations on the first differences. Conversely we may think 
of the first differences as the first summation of the second 
differences. As an example tabulate f(x) = x* for values of 
the argument from — 2 to + 4. 


Differences 
Argument | Function 
tt | 2nd | Grd | oth 
—2 —8 
+7 
—1 —1 — 6 
+1 +6 
0 Q 0 0 
+1 +6 
+1 +1 + 6 0 
+7 +6 
+2 +8 + 12 0 
+19 +6 
+ 3 + 27 + 18 
+ 37 
+ 4 + 64 
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In order to avoid ambiguity it is necessary to assign a value 
of the argument to the difference of each order as well as to 
the function. This may be done by reading along a hori- 
zontal line in the table, or along either of the two diagonals, 
We shall use all three methods and use the symbols 8, A, Vin 
designating them. They are referred to as horizontal, de- 
scending, and ascending differences. Superscripts will be used 
to indicate the order, and parentheses to indicate the argu- 
ment. Thus in the above table 


+2 =+ 8 ' 
towee WG =-t19 Wit D= 47 


pene rea +18 ve (+2) = +6 


In dealing with a small region of the table we can refer to 
a few successive values of the argument as — 1,0, + 1, +2, 
. ++: the idea being that we are not interested in the particular 
value of the argument but only in the sequence of the values 
of the function and differences. Interpolation formulae are 
usually written in this form. In this case the argument is 
written as a subscript rather than in parentheses. We list a 
few of the most common of such formulae: 


Bessel’s, 
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Stirling's, 


x (6’ 1 + 5’1/2) n* ”" 
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eGo 1) (6 1/2 +° +1/2) 


+ oe 
Lagrange’s four-point, 
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(n — 1)(n — 2) 
a Ye, 
(n + 1)(n — 1)(n — 2) 
2 = ‘ 
_(n + 1)n(n — 2) 
2 ‘ 
_ (n + 1)n(n — 1) 
6 2 
Li»j+h+in+h = i 


As special cases of Lagrange’s four- and six-point formulae we 
have for interpolation to halves: 


fin = 1/16(— fs + So + 9h — fe), 
fi 2 1/256(3f_» —= 25f-4 = 150fo + 150f; =k 25f2 =r 3fs). 


Simpson’s rule for mechanical quadrature: 


Ric = 1/3w(fo + 4/1 + fe). 
*/0 


The use of Stirling’s and Bessel’s formulae has been illus- 
trated in chapter IV. The use of formulae of the Lagrangian 
type will be discussed in section 3. The question of sys- 
tematic subtabulation will be treated in section 5. 


3. INTERPOLATION BY THE LAGRANGIAN METHOD 


The adaptation of the punched card method to formulae 
of the Lagrangian type is well illustrated by the interpolation 
to halves by the four and six point formulae. The former is 
particularly efficient. Interpolation to fourths is done by 
two successive applications of the method. We now proceed 
to describe interpolation to halves. 

Examination of the formulae of section 2 shows that in the 
6-point formula only three multiplications are required for 
each value of the function, and these are easily checked by 
addition since they involve a constant multiplier. The four- 
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point formula requires only one multiplication, for the multi- 
plication by 9 can be performed as plus ten minus one times 
the function. In the six-point formula the factor 1/256 is 
combined with the 3, 25, 150; but in the four-point the multi- 
plication by 1/16 is reserved till after the addition of the 
components. The remainder of the work consists of sorting 
the cards into proper groups and adding on the tabulator. 

Consider the case where we have three different functions 
of a similar nature tabulated for the same values of the argu- 
ment. We shall designate them as f(x), f(y), f(2), the sub- 
script indicating the argument being understood. An ex- 
ample of three such functions would be the three components of 
the force acting on a planet. For each value of the argument 
there would be a card punched in five fields: (1) argument, 
(2) f(x), 3) fo), @ f(z), (5) identification of the problem 
and of the particular step in the computation. In the inter- 
polation of perturbations for a planet we reserve columns 1-10 
for fields (1), (5); columns 11-20 for field (2); columns 21-30 
for field (3); and columns 31-48 for field (4). 

The first step after the cards have been punched and 
verified is to make one or two additional sets of cards with the 
high-speed reproducer. For the six-point formula we require 
two duplicate sets in addition to the original. For the four- 
point formula we require one set like the original except that 
f(x), f(y), f(@) have been moved one position on the card, 
which amounts to multiplication by ten. While making the 
cards on the reproducer, the new identification is applied by 
the gang punch. If the original set was identified by a zero 
in a given column we use “1”’ for the first set made on the 
reproducer and “2” for the second. In addition an “x” 
must be punched in those groups which are to be subtracted 
according to the formula, which is in both cases those identi- 
fied as “1.” In the four-point case this “‘x’’ will be used to 
add and the ‘‘no x”’ to subtract. 

For the six-point formula f(x), f(y), f(g) are now multiplied 
by 3/256 for the “0” set, 25/256 for the “1”’ set, and 150/256 
forthe ‘2’ set. This is of course done as rate card multiplica- 
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tion, and checked on the tabulator (see ch. III, sec. 2). The 
products should be so punched that when the card is reversed 
the product will occupy the same columns as the original 
number. 

The cards are now ready to be arranged in groups accord- 
ing to the proper formula. For the six-point formula we 
must make up a control group containing a card for each of 
six consecutive values of the argument, the first and the last 
being ‘‘0,’’ the second and fifth being “1,” and the third 
and fourth being ‘‘2”’ cards. In addition the tabulator must 
change the sign of the “1’’s by means of the “x” punch. 
Further the tabulator must give the argument of the result 
as that of fy plus 1/2 the tabulator interval. The groups will 
be made up with the sorter, and special “control” cards (see 
ch. III, sec. 6) will be used to give the proper control group 
and to give the proper argument of the result. 

The control cards constitute a set which serves for all 
interpolations of this type. The arguments are punched in 
the same columns as on the detail cards. In another field of 
the control card is the argument which is to go with the 
interpolated result. Thus for the six-point formula, where 
detail cards in the group have arguments 4, 6, 8, 10, 12, 14 the 
control card must sort into the deck as a three but must 
print the resultant argument as 9. One set of control cards 
with a range of one thousand in the sorting argument can carry 
in different fields the proper interpolated argument for inter- 
polation with an interval of 2, 4, etc., by the four, six, or eight 
point formula. The necessary stock of cards for this method 
is thus very small compared to that of any other method. 

The sorting into control groups is done as follows for the 
six-point. All the cards ‘‘0,” “1,” “2” are sorted together 
according to the argument, and separated into packs accord- 
ing to the period six (see ch. III, sec. 4). They are then 
further subdivided according to the “0,” “1,” “2.” There 
are thus eighteen packs of cards which are placed on a table 
in three columns of six rows each. The three columns 
contain, respectively, the ‘‘0,” “1,” and ‘‘2”’ cards, and the 


= aia beens nny ome Se od 
Ne aa 


ee 


Speer - 
SS leer ey ha ee eee als 


iss 


coh ca 


ri 


an 


i 
Hy 
a 
| 
.| 
3 
} 
) 
i 
t 


50. Ss Se : ‘PuNCHED Carp MrrTHops 


‘ix rows contain cards whose arguments, after division by 6, 
ave ceadinalll 0, 1, 2, 3,4 5. If now we take up from the 
table all those whose position is indicated by “a” in figure 13A 
we have all the cards necessary for obtaining every sixth 
interpolated value in the table. These cards are placed in the 
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sorter together with the set of control cards, and sorted in 
order of the argument. This set is now ready for tabulation. 
When this set has been tabulated and the results recorded 
with the summary punch, we proceed to do the same for the 
cards whose positions are indicated by “b”’; and finally by ‘‘c.” 
It is noted that the scheme of figure 13A gives half the 
required groups. The others must be gotten by replacing the 
cards on the table and collecting according to figure 13B. 
The above sorting may be done at one time for any number 
of sets of functions having the same arguments, a final sorting 
being required to separate them. 

~ For the four-point formula are required four consecutive 
values of the function “0,” and two values of “1.” After 
the summary cards have been obtained the interpolated values 
must be multiplied by 1/16: 


4. Tue ForMATION OF DIFFERENCES 


The formation of successive differences of a function is 
the first step in most methods of interpolation and numerical 
integration. We give here twoalternative methods of forming 
these differences with a tabulator. In the first method the 
differences of all orders are derived for a given value of the 
argument before proceeding to the next, while in the second 
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the differences of the first order for all the arguments are 
derived before proceeding to the second. The first method 
is necessary in numerical integration, and may be desirable 
where high-order differences are involved, especially if the 
table contains errors which are to be found from the differ- 
ences. The second method is more efficient in that it does 
not require the full attention of the operator. 

In either method all the counters of the tabulator and a 
large part of the card should be used at each step. If there 
are several functions of a similar nature for each value of the 
argument, they should be taken together. For instance, in a 
great deal of astronomical work we use the same card for f(x), 
f(y), f(z). If there is only one function and the table is long, 
it is possible to break the table into three or four parts so 
that the values of the function for arguments separated by 
100 or 1000 may be considered together. The first counter 
of the tabulator is used for the argument, the order of the 
difference, and the identification of the particular set of func- 
tions being differenced. The second, third and fourth coun- 
ters carry the corresponding values of the differences of the 
three functions. In some cases it might be desirable to use 
the four counters for four functions and to apply the argument, 
order of the difference, and identification by means of the 
high-speed reproducer. 

In the first method we wish to find all the ascending 
differences V‘ for a given argument before proceeding to the 
next. This is accomplished by finding first the difference of 
the highest order and then obtaining the others by successive 
additions. By definition we have for any diagonal in the 
table: 


CE. eC le Bua er ne A LS 


Adding, we find 
Vom! = fo — (fa + Vil + V1” + ++ + V8") 
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Thus at each step of the calculation we subtract all the num- 
bers along the diagonal from the next value of the function. 
We then add to this result the same previous numbers one 
after the other (progressive totals). This method is obviously 
suitable for the punched card technique. It is only necessary 
to work out a simple routine which will make the proper ad- 
ditions and subtractions, and to punch summary cards and to 
clear the counters at the proper times. The manipulations 
of the cards and control switches must be made simple. 

The tabulator is wired so that when the cards are tabu- 
lated in the direct position the numbers are added (alge- 
braically), and in the reverse they are subtracted. For this 
purpose an “x” must be punched in one column of all the 
cards. The dial switches are set for “total, non-reset,’’ and 
the counters are cleared by tabulating the last summary 
card reversed. Minor control is wired to the column which 
indicates the order of the difference, and the corresponding 
column on the reversed card must have a zero in each card. 
Three colored cards, the first punched zero in the control 
position, and the other two punched with the digit of the high- 
est order difference, facilitate the routine. 

The operator’s routine is: 


~ (1) Turn switch “S. P. Minor” to “off” position. 


Q) Tatulate: fo, fa, Ya, V-', ---, V-" reversed; fo, 
three control cards direct. 
(3) Place in tabulator hopper: V_:", Vi", «++, V-i’, fu. 


(4) Turn switch “S. P. Minor” to “on” position, and 
start tabulator. 


In the above routine we have started with all the values of 
fo punched on cards, and have in addition formed a duplicate 
card for each value to simplify the clearing of the counters. 

In the process of numerical integration (see ch. VIII) we 
form five or six differences and two summations for f(x), 
f(y), f(@) by an adaptation of the above method. For that 
work a special circuit in the calculation control switch pro- 
vides for the clearing of the counter, thus eliminating the 
additional card fo. 
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In the second method we wish to subtract (algebraically) 
the numbers on each card from those on the following. This 
requires that the cards be tabulated in pairs, and that the 
first card of each pair be subtracted and the second added. 

The necessary wirings for both these operations are de- 
scribed in appendix B. Section D of figure 21 refers to the 
tabulation of the cards in pairs. Section A gives the wiring 
for subtracting numbers on one set of cards from those on 
another. The subtraction requires an ‘‘x’’ in alternate cards. 

Tabulation of the function cards with the above wiring 
gives the first difference for every other value of the argument. 
When this operation is completed we remove the first card of 
the pack and run them through again to get the remaining 
arguments. For the second run we must change the wiring 
of the ‘‘«” so that in this case also we subtract the first 
card of each pair from the second. 

For inspection of the differences it is necessary to sort 
the two sets of differences together and list them. 

The ‘‘x”’ on the function cards may be put on by the key 
punch operator or with the gang punch. The “x” on the 
first differences, needed to get the second differences, may be 
applied during the summary punch operation. 


5. THE SUBDIVISION OF TABLES 


In the construction of tables where each value of the 
function involves long computation, the function is computed 
for a few key values and the ones between are inserted by 
interpolation. We have discussed the Lagrangian method 
of interpolation to halves and quarters. Where more than 
three values of the argument are to be inserted the difference 
method has a great advantage. We have seen in chapter IV 
how tables may be easily built by summation, once the differ- 
ences are know. The problem is then to find a method for 
determining the differences in the subdivided table from those 
of the original. To be suitable for the machine the method 
must be continuous, i.e., it is not desirable to make any 
numerical adjustment by “‘inspection.’’ On the other hand 
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a machine method could be efficient even if more arithmetic, 
within reason, were involved than for a hand method. 

Two methods have been suggested for use with punched 
card machines, and I shall refer to them as the ‘end figure” 
method, and the “superfluous difference’ method. Both 
methods are important contributions, but from the nature of 
the problem, it seems that a still better method could be 
devised. 

In the end figure method the end figure of each inter- 
polated value of the function is computed. These end figures 
are differenced to an order high enough to give smooth dif- 
ferences, say second or third. Where the differences run 
smoothly the entire value of the difference can be filled in 
when the end figure is known. With these values of the 
highest order differences, those of the lower orders are built 
up by summation. The end figures of the interpolated values 
of the function are formed on the tabulator. Cards are se- 
lected from the files depending upon the various order dif- 
ferences of the original table, and these are added on the 
tabulator. (See reference (3), p. 43.) 

The end figure method has the advantage of giving the 
interpolated functions and differences just as they are to be 
printed in the final table. It has the two disadvantages: 
(1) a collection of thousands of special cards is required for 
the interpolation; (2) a large amount of clerical work is neces- 
sary in selecting the proper cards from the files, and in 
differencing the end figures. This clerical work could of 
course be adapted to the machines if sufficient interpolation 
were contemplated to warrant it. 

The superfluous difference method builds the subdivided 
table from differences of one order higher than that required 
in the original table! In addition to the high order of the 
differences, extra decimal places are required for the summa- 
tion. Another disadvantage is that, having used extra‘ deci- 
mals in the differences, they will not be the values of the 
differences which should be printed in the table. After the 
functions have been properly rounded, the end figures of the 
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differences must be adjusted according to the rounded end 
figures of the functions. From reference (5) at the beginning 
of the chapter is given below a table showing the differences in 
the subdivided table in terms of those in the original. It 
must be remembered that in building tables from any given 
difference by this method, the accuracy of the interpolation is 
only that of one order less. 


Differences 
Third Fourth Fifth Sixth Seventh 
+.055 — 0165 — 0086625 +.00329175 + 0017007375 
—.045 +.034 +.0257125 — 0134145 — 0085259625 
0 — 0165 — .0247875 + .0202555 +.016953475 
0 0 +.0078375 — 0134145 — .016716525 
0 0 0 + 00329175 + .0081802875 
0 0 0 0 — 0015910125 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 i 0 


For comparison of the end figure method with the super- 
fluous differences method let us assume that interpolation 
must include fifth differences, the sixth being negligible. In 
the superfluous difference method we must build with sixth 
differences using eight extra decimal places. In the end 
figure method the third differences in the subdivided tables 
are of course only about one thousandth of those in the 
original, and no extra decimals are used in the summation. 
With ten-digit counters in the tabulator and six-digit func- 
tions, the superfluous difference method would work well for 
fourth. differences. The alphabetic machine has greater 
counter capacity and could be used for this purpose. A 
numerical tabulator with special carry-over from one counter 
to another would be useful. 

In all these difference methods one is impressed by the 
number of cards which must be punched by the machines and 
used only once or twice. In England there is on the market a 
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machine which will transfer from one counter to another, 
This “rolling total"’ machine would be much more efficient for 
forming differences and summations since it would not be 
necessary to punch cards. This observation does not apply 
to the process of numerical integration for then the tabulator 
is used for other operations in the meantime and the numbers 
could not remain in the counters. 

The newer models of the multiplying punch offer possi- 
bilities in the building of tables by summation. In these 
machines it is possible to transfer from one counter to another 
and to punch progressive totals. Withslight modification one 
of these multipliers might become a remarkable ‘ Difference 


Engine.” 


CHAPTER VI 


NUMERICAL HARMONIC ANALYSIS 
AND SYNTHESIS 


1. INTRODUCTION 


A great many functions occurring in science may be repre- 
sented with sufficient accuracy by harmonic series containing 
a limited number of terms: 


(1) F(x) = a) + a, cosx + a,cos2x +---+a, cos nx, 
+ b,sinx + by sin 2x +---+6,_,sin (n — 1)x. 


The problem of analysis is the determination of the a’s and b’s 
from special values of F(x). The problem of synthesis is the 
determination of special values of F(x) from a given series. 

The usefulness of the method depends to a great extent 
upon the ease with which the analyses and syntheses can be 
carried out. Both are done with great efficiency by the 
punched card method where the series have many terms or 
where many functions are to be analyzed or synthesized at 
the same time. 

The numerical work in either case consists of the solution 
of linear algebraic equations, and the special values of x are 
chosen in such a manner as to simplify the numerical work. 
If we choose special values of x as follows: 


x; = th, h=rn/n 
we shall have the well-known formulae: 
n—l 


(2) | = a, cos kh + ¥ by sin ikh 
k=l 


k=0 


(¢ = 0, 1, 2, --+, 2n —4), 
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nea; = "SF; cosijh (7 =9,1,-°*, %), h = x/6 ip i 
va d=sinh=¥% €= cosh = Yv3 i 
6 =2 for j=O,n 1 
= 1 fee ge 1,2, --*,e — 1, Sin it 
2n—1 P; oe _ SS ee ee ae a4 
a F; sin tjh G= Beds 2s 1). j at 
For a great many problems it is more convenient to sepa- : Bear at 
rate the sines from the cosines and consider each separately, ) 0 0 0 0 0 0 a if 
Let I O +d +e +1 +e +d 0)2 : 
ea CS, 2 Dae. — AN’ 0 -e -e O0Oj2 v 
where 3 0 +1 aes eee oe oe i 
+B 
C(x) = a + a cos x + +++ + Gn COS MX, 4 ie | eee ee } 
S(x) = b sinx +---+ d,.sIin (n — 1)x. ; ee ee eS di 
Then for values of x chosen in pairs, a and 27 — a, we have 6 0 0 0 0 0 0 0 | i 
| 16s 
F(a) + FQx — a) = 2C(a), 6 6 6 6 6 4k 
Fla) — F(2Qx — a) = 2S(a), i 
and the formulae of analysis and synthesis become: a 
(3) C= za cos ikh a> 01, --,”), | a 5 ; ; : : : : 
al : Se 
S; = Xb, sin kh G=1,2,---,n — 1), ¢) aes | ae! | Sa SS, so] a] Al 1 
k=1 
: i +1 +e +d 0. = d- =e = 2 
nea; = SeiC; cos ijh G=0,1,--:,7), 2 Sal ao =a = 1 Sd See ee 
- 3 + 1 QO -1 0 +1 QO -1 2 
n—l 
nb; = > «,S; sin ih 1,2, ---,n — |), 4 ne | Lee ag et et | =e a | 2 
i=1 
5 —e = 0 -d —1|2 
qed fers. 4 = 0, 7 +1 eo 48 oe 
= ioe se = 1,2, «>, n — 1. 6 eee Din tee PG Ela Ee Ee ee 
The analysis or synthesis is easily executed with the help ne; 12 6 6 6 6 6 re ai 


of a table of sin (ijh) and cos (ijh). We give here such a 
table for n = 6. 
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Each computation consists of multiplying given numbers 
by the numbers from these tables and adding the results, 
A great deal of the work consists of adding the same numbers 
+n different combinations, which is done most efficiently by 


punched cards. 
2. MacuINE METHODS FOR SMALL SERIES 


The machine technique is best described by means of an 
example. Consider the analysis of a cosine series when = 6 
and where each value of C; is given to seven significant 
figures. Suppose there are 99 such sets of functions to be 
analyzed at the same time. For each set are given C(0°), 
€(30°), C(60°), C(90°), C(120°), C(150°), C(180°) or C; where 
i=0,1,---, 6. 

Since our tabulator has four counters it will be possible 
to add corresponding quantities for three different functions 
in three of the counters and to use the fourth for identification 


of the series and the operation. Nineteen cards will be re- 


quired for each three sets of analyses. 

The given functions are punched on the cards as follows: 
Columns 1-2 contain the number of the series being analysed, 
i.e. 1, 4,7, ---, 97. Column three contains z = 0, 1, ---, 6. 
Columns 5-11 contain C; for the series whose number is 
punched in columns 1-2, i.e. numbers 1, 4,7, ---, 97. Column 
12 is punched with an “x” if C; is negative, and with a ‘“‘y” 
if positive. C; for numbers 2, 5, 11, ---, 98 are similarly 
punched in columns 14-21, and for numbers 3, 6, 12, ---, 99 
in columns 23-30. A zero is punched in columns 4, 13, 22. 
When the cards have been punched and checked they are 
passed through the reproducer to give the necessary dupli- 
cates and to punch quantities needed for the tabulation. It 
is necessary to make duplicate sets fori = 1 — 5. To iden- 
tify the duplicates a “‘y”’ is punched in column three. All 
cards with 7 = 1 and 5 (including duplicates) are now sorted 
on columns 1, 2, and the same numbers gang punched in 
columns 80, 79, Finally we gang punch on the cards the 
following pattern: 


Gi ani 


aes rr se 


SLA Pe Ee ST 


NuMERICAL Harmonic ANALYSIS AND SyntHeEsIs 61 


Column on Card 


Value of 

t 
31 | 32 | 33 | 34} 35 | 36 | 37 44 |45 46 47 | 48 49 | 50 

0 Gebers law (5-16 
1 0 2 Xs Xs Xs i 
ly 0 Xe | 
2 0 1 ta Xs x, 5 6 
2y 0 Xs 6 
3 0 X 4 xs 
3y 0 Xa 4 x | 
4 O |} m | | 3 | % | m1 6 
4y 0 3 6 | 
5 0 2 Xs x6 5 x 
Sy 0 Xs 
6 Celene Xe i 4 | x16 


| 


This pattern insures three things in the later tabulation: (a) 
that the proper cards are-included in each group, (b) that 
they are tabulated in the required direct or reverse position, 
and (c) that each has the proper sign. 

The numbers on cards 1 and 5 are now multiplied by 2e 
as rate card multiplication, and checked on the tabulator 
(see ch. III, sec. 2). The products are punched so that on 
reversing the card from left to right the product will occupy 
the same position as the multiplicand. No signs are indicated 
in the products. Three wirings of the multiplier are necessary 
because three functions occur on each card. 

The tabulator is wired to control on columns 1, 2, 31. 
Counter (1) is wired as follows: positions 10, 9, 8 to columns 
1, 2, 31 of card through group indication; position 7 to card 
count. Counters (2), (3), (4) are wired to columns 4-12, 
13-21, 22-30 respectively. 

Counter (1) always adds. Counters (2), (3), (4) add and 
subtract from ‘‘x” in columns 12 or 69, 21 or 60, 30 or 51 
when there is no x in 31, and subtract and add when there is 
an x in 31. This is accomplished by the method described 
in chapter II, section 3. 

The summary punch is wired to punch positions 10, 9, 8 
from counter (1) in columns 1, 2, 3 on the card. The quan- 
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tities from counters (2), (3), (4) are punched in columns 4-12, 
13-21, 22-30 with “x” in columns 12, 21, 30 for negative num- 
bers. This summary card therefore resembles the original 
cards. : 

All the cards are now sorted in order according to columns 
1,2. Tabulation with summary punch on minor control gives 
12a» for each function. The results are checked by tabulation. 

The cards are now sorted oni. The cards fort = 0, ly, 4, 
6 are put in the sorter in the usual manner, but those with 
i = 1, 5 are reversed. All are now sorted on columns 1, 2. 
The plug wires on the tabulator going to column 31 of the 
adding and control brushes are moved to column 32 and every- 
thing is ready to give 6a. 

The above process is continued according to the schedule 
until i = 6 has been completed. 

We have thus obtained 12a, 6a:, 6a2, ---, 125. Rate 
card multiplication by 1/12, 1/6 gives do, a, «++, ds. These 
results may be listed or used for further computation. If the 
results are to be synthesized they may be run reversed through 
the reproducer to obtain a set of cards of the same form as the 
original. 

The modifications necessary for the synthesis or the analy- 
sis of a sine series are obvious and need no further description. 

The checks outlined above are very rigorous, and the 
general appearance of the results is usually a good test of the 
work. However, in very important cases the results in analy- 
sis may be checked by synthesis and vice versa. 

The time required for this example should be: 


Punching and checking............... 3 hrs. 
Reproducing and gang punching....... Y hr. 
Re ee a css ..- 2 hrs. 
Tabulating and summary punching..... 2 hrs. 


3. Macuine Metuops ror LARGE SERIES 


In the above example the machines were used advan- 
tageously on series of only a few terms in those cases where 
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many analyses or syntheses were to be done at the same time. 
For series of many terms it may pay to use the punched card 
method for only a few series at a time. 

Consider for example the cosine series with thirty har- 
monics where numerical values of the coefficients are given 
and it is required to compute C(x) for the special values of 
x = nh, h = 6°,n = 0,1, 2, ---, 30. Four figure accuracy 
will be assumed, though this can be anything up to eight. 

For this method a set of 228 prepunched cards are used 
for each three series to be synthesized. They contain the 
required cosines and also 30 columns with numbers for sorting 
and for giving the proper algebraic sign of the products. The 
a’s are punched by hand onto other cards. Each “‘a”’ card 
will contain corresponding a's for three different series. Nega- 
tive values are indicated by an ‘‘x”’ punch. The “a” cards 
and the prepunched cards are sorted together and passed 
through the multiplier to give the product of each a by the 
required cosine; the ‘“‘a’’ cards are the rate cards and the 
prepunched cards are the detail cards. The multiplication 
is checked as it is completed by means of the sum of the 
products (without regard to sign). The cards must go 
through the multiplier three times, once for each of the three 
series. In addition to the above check on the absolute values 
of the products, the sign values are checked by one run on 
the tabulator. The cards are now sorted according to the 
value of ‘‘x’’ and added on the tabulator. If the tabulator 
has a summary punch the results may be checked by addition. 

We have set aside 30 columns in the prepunched cards for 
sorting, control, and indication of the signs of the products. 
These 30 columns correspond to the desired values of “7.” 
To add the products for a given value of “‘n”’ the cards are 
sorted on that column and all the cards which are not to be 
added are rejected by the sorter. An ‘‘x’’ in this column is 
used to operate an “x-distributor’’ on the tabulator in such a 
manner that terms with negative a’s will be added, and those 
with positive ones will be subtracted. 
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The time required to synthesize 48 series will be: 


Reproduce the prepunched cards (16 


Te Ie So Sao 1 hour 
Punch and check the oat Se ee 3 hours 
Muley xe cere tc ttt 15 hours 


Sort, tabulate, and summary punch.. 10 hours. 


The methods described above in this and the previous 
section require very little in the way of prepunched cards, 
and hence leave great freedom in the choice of the particular 
scheme to be adopted. Should occasion arise where one 
particular scheme is to be used day in and day out, the stock 
of prepunched cards could be efficiently increased. For ex- 
ample Dr. Linus Pauling of the California Institute of Tech- 
nology in his work on crystal structure has synthesized 
thousands of series of thirty terms each. It has proved 
efficient for his work to prepare a set of about 65,000 cards 
containing all the possible combinations of coefficients. He 
can then select the proper cards by hand and add the com- 
ponents on the tabulator. 


CHAPTER VII 


THE MULTIPLICATION OF SERIES 
1. THe PRoBLEM 


In many physical problems we have two functions each 
expressed as a harmonic series with a limited number of terms 
and we require the harmonic series which represents the 
product of the two functions. The formation of these product 
series is efficiently done by the punched card method when 
the series are long or when many product series can be formed 
at one time. The application of the method here described 
to power series is obvious. 

We shall discuss the case where the argument of each 
term is a linear function of the independent variable involving 
several fundamental frequencies. This includes of course the 
ordinary Fourier series. Three kinds of series occur: sines, 
cosines, and exponentials. 


Let 
(1) S(x) = 2b; sin ai, 
C(x) = Ya; cos aj, 
E(x) = Zc,e'%, 
where 
a; = (jp + kg + ir + ---)x, 
j, eb, +> =O, Lb 2 ss, 
~, 9, 7, +++ = given numbers. 


The product of any one of these by itself or by any one of the 
others may be required. In forming the product of two series 
it is necessary to consider each term of the first in connection 
with each term of the second. The two coefficients must be 
multiplied and the harmonic functions combined according to 


65 


66 Puncuep Carp METHODS THe MuLrtiPiicatTion oF Serres 67 


the trigonometric formulae: Since the argument has more than one component it would be 


a) danas g wen(e— 2) cate +) Bee cartons pen oe i pati 
eg ean te 9) + in (x + y), of si is = we decide arbitrarily on the proper 
Fey eee y) + cos (« +9), order o the components of the argument and treat it as 

orev = ert, positive or negative according to this scheme. Thus we may 

Sg aaa ees, | consider an argument as positive if m is positive and negative 

a eee oe. if m is negative. If m equals zero, | decides, and then & and 

- : finally 7. This leaves only the constant term. The coeffi- 

Each pair of terms in the product of two trigonometric cients are punched as true figures with an ‘‘x” to indicate 
series gives a term with the sum of the arguments and one negative numbers. The integers 7, k, 1, m however are made 
with the difference. Exponential series require only the sum, positive by the addition of say 50. Thus the argument 

but since exponential series are usually used for the computa- . (p — 2qg + s) would be indicated on the card as 51485051. 

tion of functions where both the function and its conjugate When the cards have been punched and verified they should 

are required, the difference of the arguments is needed also be listed in the order of the argument for proofreading. 

for that case: i The cards of each series are now sorted according to the 

The multiplication of one harmonic series by another thus numerical value of the coefficients with 1 « 107° first and 
consists of forming all the sensible products of the coefficients ; 1.0000000000 last. It may seem at first thought that this 
and the corresponding sums and differences of the arguments, E would require ten complete sortings, but this is not the case. 


and collecting these terms according to the resultant argu- PB The chance that a number is zero in the units and tens posi- 
ments. 
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tion and not in the remaining positions is only one in one 
The machine technique will vary with the particular prob- hundred. Hence after the second sorting large groups may 


lem, depending upon the number of terms in a series, the be eliminated from further sorting by looking through the 
number of fundamental periods in the arguments, and other : holes. After the third sorting the total number may be still 
considerations. We shall describe the procedure for two prob- : further reduced. The final sortings involve only a few cards. 
lems of widely different nature. One or the other of these We now print with the tabulator two lists for each series. 
two methods should be adaptable in most cases. The first has four columns: (@) consecutive numbers from 1 
to 500, (6) the argument, (c) the coefficient with sign, (d) the 
coefficient without sign. The second list is just like the first 

Let us consider the case where we have S(x) and C(x) except that progressive totals are printed in columns 8, ¢, d. 
each with five hundred terms, four fundamental periods in It is now necessary to set a rejection limit for the products, 
the argument, and ten significant figures in the largest co- b say 1 X 10-".. In other words all products smaller than 
efficient. Suppose for simplicity the largest coefficient is 1.0. ) 1 x 10-" will be rejected. We must examine each pair of 
Assume further that extreme care is to be exercised to prevent coefficients to see if the product is within the rejection limit. 
error. In such a series the integers j, k, 1, m would not Consider the largest term of C(x); it will give a sensible 
exceed + 10. product for each term of S(x). The next largest term in C(x) 

Each term of each series is punched on a card showing may permit the rejection of all terms in S(x) smaller than 
identification of the series, the argument, and the coefficient. say 2 X 107°. This rejection limit is computed for all terms 
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in C(x) greater than ¥10 X 10-". It is more convenient t 
examine the remaining terms in the reverse order, i.e, pi 
pute the rejection limit for the largest term in S(x) and so 0 
until we come to the term in S(x) whose rejection limit ‘ 
Vi0 X 10-*, Thus all the sensible products may be foy 
by considering about 100 terms. ee 

_ With the help of the computed rejection limit we form all 
the required pairs of terms. A new card is made with th 
high-speed reproducer for each pair of terms. These a : 
made as follows: the first card of C(x) must be paired wi 
each card of S(x), the second with all greater than the rejec 

tion limit, etc. The set S(x) is placed in the feed A of a 
reproducer; and new cards, preceded by the proper card fro : 
C(x), are put in feed B. Thus by a combination of a 
punching and reproducing we obtain cards each of hick ue 
two arguments and one coefficient. The coefficient which : 
common to the whole group need not be punched on ie 
cards but may be punched on a special red rate card to go 
with the group. By a continuation of this process we obigi 

about 100 groups of cards each with its rate card. Since the 

capacity of the multiplier is eight digits, those groups which 
have 9 and 10 digit multipliers and eight digit multiplicand 
are done in two groups. The few terms with large multiplier 
multiplicand are done by hand. In a series with ee 

‘ ~ convergence there may be about 20,000 cards altogether, 

: reproducer is of course used to check the punching. 

tee to the check of the reproducer itself we may 

: . an a tional check from the tabulator. Each group 

ae is tabulated to give the sum of the arguments, the 

prin the sum of the coefficients, and the sum of the 
re Aas of the coefficients. These totals are now com- 

Be, a : ie “a eae by subtracting the proper 

Hes with the formation a ae oS 

aS — ee = eee performed for each group and 

Ley ao iplication. Extra digits are accu- 

€ summary products counter as described i 
unter as described in 
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chapter III, section 2. The check sum from the summary 
products counter is compared at once. This sum checks the 
multiplication but not the signs. Additional checks are ob- 
tained during a later tabulation. 

The value of the sum and difference of the two arguments 
on each card is evaluated by means of the sorter and high 
speed reproducer. For those cases where the value of j, B, / 
or m which is to be added and subtracted is zero, the desired 
result should be applied at the time the duplicate cards are 
made; it being only necessary to punch the same number for 
the argument and for the sum and difference of the two 
arguments. The largest groups of course have mostly zeros 
in the argument, and hence most of the work is done at the 
time the duplicates are prepared. For the others we have a 
set of master cards giving the sum and difference of any pair 
of integers from minus ten to plus ten. These cards are 
sorted together with the detail cards for 7 of the two argu- 
ments, and run together through the gang punch. The wires 
are then moved on the reproducer and the cards are re-sorted 
to give the corresponding results for k and then / and m. 

When all the sums and differences of the arguments have 
been formed the cards are sorted to test for the proper sign 
of the product, and are then sorted into the original groups. 
These groups are now tabulated to give the sum of the dif- 
ferences of the arguments, the sum of the sums of the argu- 
ments, the sum of the products of the coefficients, and the 
sum of the absolute values of the products. The sums and 
differences of the arguments are checked exactly. The other 
two sums give fairly accurate checks on the punching of the 
products. 

The cards are now sorted in order of the sum of the argu- 
ments and tabulated with control on the sum of the arguments. 
The results are recorded on the printed record and on the 

The resultant argument of course comes 
through the group indicator. The checking is done by tabu- 
lation of the old cards and the summary cards. We have 
outlined checks for such work in chapter III, section 2. 


summary cards. 
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Here every precaution should be taken with the checking, 
such as reversing the cards in the tabulator, and interchanging 


the counters. 
The cards are now sorted according to the difference of 


the arguments, and then separated into two groups according 
to whether the cards were made with C(x) or S(x) as common 
multiplier. In other words the difference of the arguments 
may be (x — y) or (vy — x) of equation (2). The summary 
cards are obtained for the two groups separately, and by 
means of the switch on the summary punch the signs of the 
coefficients are changed for the set (y — x). 

We now have three sets of summary cards with arguments 
(x+y), (x — 9), (9 —x). They must be combined, and 
at the same time those having negative arguments may be 
combined with those having positive. Since the resultant 
terms are sines, the sign of the coefficient must be changed when 
the sign of the argument is changed. We sort (x — y) and 
(y — x) according to the sign of the argument. For negative 
arguments we punch the complement of each component with 
the aid of the sorter and reproducer. These new arguments 
are punched so as to occupy the same position when the card 
is reversed. The cards (x + y), (x — y), (vy — x) are now 
combined, those with positive argument direct and those with 
negative argument reversed, and sorted together according 
to the argument. The cards are tabulated and the results 
summary punched to give the resultant series 2 C(x) S(x) ex- 
pressed as a sine series. 

In forming the square of a series only one stack of cards 
is involved, and fewer products are necessary. Each card 
need be taken only with itself and those following. This 
gives all the distinct products but not all the products of the 
required series. Since there are many more cross products 
than square terms the easiest thing to do is to divide all the 
squares by two. This gives the square of the original series 
for sines or cosines, and half the square for exponentials. 

The following identification has been used to identify the 
various cards used in forming a product series. Two columns 
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are required for identification. The first contains: 


0 = first original series, 
1 = second original series, 
= red rate cards made from first series and the correspond- 
ing detail cards made from the second series, 
blue rate cards made from second series and the corre- 
sponding detail cards made from the first series, 
= summary cards (x — y), (y — x), 
= summary cards (x + y), 
= product series. 


second contains a zero except: 


(1) Groups 3, 4 for the large digits of those groups whose 
multipliers must be separated into two parts. 
(2) Group 5 for (y — x). 


In these cases a “‘1” is punched instead. In order to secure 
uniformity those product terms which must be done by hand 
are punched entirely on the (03), (04) cards, and no (13) or 
(14) cards are used for these terms. Also, as soon as all the 
products have been punched and checked the square cards 
are removed and new ones made on green cards with half 
the coefficients. If the sum of the products is formed first 
for (x — y) these terms all come in one group with arguments 


50505050. 
3. SMALL SERIES 


The case outlined in section 2 is an example of the largest 
series that are apt to be encountered. One such series is 
sufficient to justify the use of the punched card method. For 
small series it would be necessary to do a considerable number 
at one time, unless the previous and following computations 
are also to be done on the punched cards. 

Consider the case where we have one hundred pairs of 
series to be multiplied. Each series has one fundamental 
frequency and nineteen values of ¢ from plus 9 to minus 9. 
The coefficients have eight digits, and products less than 
1 X 10-* are to be neglected. 
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The cards are punched so that each card contains the 
argument, the number of the series, and the coefficient. 
Usually in such a case as this some of the series will each be 
used to multiply several series, but only one set of cards 
will be punched by hand. 

When the cards have been punched they should be listed 
with the tabulator and totals obtained for the sum of the 
arguments, the sum of the coefficients, and the sum of the 
absolute values of the coefficients. 

Duplicate sets should be made with the high-speed repro- 
ducer so that for each product series required there will be a 
multiplier series and a multiplicand series, each carrying a 
number to indicate the number of the product series. Colored 
cards should be used for the multiplier series. The multi- 
plicand cards should have an “x” in a given column. In 
case only a few duplicate cards are required this step could 
have been avoided by modifying the initial hand punching. 
All the sets of cards should now be listed as for extensive 
series, and the totals punched on summary cards. 

For these small series the constant multiplier groups are 
so small that the hand method of rejecting the negligible 
products would not be efficient. It must be done by succes- 
sive sortings, on all the series at one time. 

The colored multiplier cards for all the series are now 
sorted according to the numerical value of the coefficient, and 
the corresponding “‘rejection limit’’ for terms of the multi- 
plicand is applied with the gang punch. This latter quantity 
need be accurate to only one or two significant figures, but 
should be rounded so that additional terms will be included 
in each group rather than that significant terms will be 
omitted. It should be applied so that on reversing the cards 
it will occupy the same field as the original coefficient. The 
argument and number of the series should also be transferred 
to the other end of the card so that when the card is reversed 
they will occupy the same columns as before. 

The multiplier and multiplicand cards are now sorted 
together according to the numerical values of the multi- 
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plicand cards. The multiplier cards are placed in the sorter 
in the reverse position followed by the multiplicand cards 
in the direct position. When the cards are all in order accord- 
ing to the coefficients, they are then sorted into 100 groups 
according to the number of the series. In each of these series 
the multiplier cards (red) are followed by the necessary multi- 
plicand cards, and the corners of the multiplier cards stick 
out so the top group in each case may be easily lifted from 
the pile. These top groups are collected and the red cards 
s2parated from the manila ones by sorting on the column with 
the ‘‘x.”” The manila cards are now put through the repro- 
ducer to form duplicates, and the duplicates are then placed 
with the red ones. The cards which were used to form the 
duplicates are replaced in the reproducer and to them are 
added all those groups obtained by lifting the second colored 
cards from the original series. These second colored cards 
are attached to the second set of duplicates as before. This 
procedure is continued until all the necessary duplicates 
have been formed. Finally there will be nineteen packs of 
cards containing red and manila cards. These are all placed 
in the sorter, in each case red cards followed by the corre- 
sponding manila ones. They are all sorted in order of the 
numbers of the series. By means of the gang punch we now 
transfer the argument and the coefficient of each red card 
to the following manila ones. The “x” punch is used to 
prevent the transfer to the red cards. The red cards are then 
separated from the manila ones by means of the sorter. We 
now have all the cards of a given series ready for the multiplier. 
There should be perhaps forty or fifty in each series. 

The multiplication for each series is performed, and the 
products punched with nine digits. The sum of the products 
in the summary products counter may be checked as each 
group is completed by comparing it with the product of the 
sum of the absolute values of the multipliers by the sum of 
the absolute values of the multiplicands for the series. 

With the tabulator the sum of the punched products may 
be formed. Each of these should be compared with the 
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product of the sum of the multipliers by the sum of the 
multiplicands. If the check products were formed on punched 
cards the comparison can be made with the tabulator. 

All the cards are now sorted on the two arguments and 
the sum and difference applied with the gang punch. 

From here on the procedure is the same as for large series. 

The machine time for 100 series is as follows: 


Original punching and verifying 
Listing and .summary punching check 

a 

4 day 
Sorting and gang punching duplicates. . 1 day 
Multiplying 1 day 
Sorting, tabulating, and summary 

punching 


In this example the checks would have been exact had 
we included a card for each product regardless of size, but 
this would have required about forty thousand cards instead 
of four or five thousand, and the multiplying time would 
have been over a week instead of one day. 


CHAPTER VIII 


THE NUMERICAL SOLUTION OF DIFFERENTIAL 
EQUATIONS 


The method described here for the numerical solution of 
differential equations will work for the vast majority of cases 
arising in physical science. It includes all simultaneous sys- 
tems of ordinary differential equations where the accelerations 
are functions of the time and of the coordinates and their 
derivatives. 


Let x, y, 2 = the dependent variables, 


t = the independent variable, 
dx ae ; 
x= WP’ and similarly for Y, Z. 
The differential equations are 
(1) X = X(x, 9, a5 4’, 9’, 25 b) 


and similarly for Y, Z. Higher derivatives may also be 
included. The solution is given as a table of values of x, y, 2 
for equidistant intervals w of ¢. It is given by the well-known 
expression 


Q) «="X + 1/12 X — 1/240 X" + 1/1951 XIV +--+, 


and similarly for y,z.. Here ’"X, X, X”, X'Y, --- are numbers 


- taken from an ordinary difference and summation table. “X 


is a column of numbers of which X is the second difference, 
X”’ is the second difference of X,and soon. The six arbitrary 
constants of the solution may be taken as the values of “X, 
"VY "'Z for any two successive tabular values of ¢. 

All the numbers in equation (2) refer to a given value of 
t, i.e., the values of X”, X, "X, +--+ lieona horizontal line in 
the difference table. In the notation of chapter V they are 
gx, O° Xe, OF UX, 
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The solution is built up step by step as follows. As soon 
as a new value of X is computed another diagonal line in the 
difference and summation table may be found. This gives 
the next value of '’X, required by equation (2), but not the 
required values of X, X", X1Y, +++: These latter must be 
extrapolated in the table with sufficient accuracy for this 
purpose. X must be extrapolated one tabular interval, X” 
two, X!Y three, etc. If the interval of tabulation is small 
enough this extrapolation may always be done correctly. The 
difference between the extrapolated value and that computed 
later is evident from the difference table. 

The values of x, y, s just derived from equation (2) are now 
inserted in equation (1) to compute the next value of X. The 
values of x’, y’, ’ for the same value of the argument are also 
required, and are computed from 


(3) x’ = (‘'X) —1/12(X") +---, 


where ('X) and (X") are the means of the numbers lying just 
below and above the horizontal line of "X, X, ---. 

The mechanical procedure is as follows. Special tables of 
functions required for the computation of X (x, y, 2; x’, y’, 2’; 0) 
are prepared in advance in the form of files of punched cards, 
each card containing argument, function, and difference. All 
the multiplication and interpolation of tabular functions re- 
quired for equations (1), (2), (3) are done on the multiplying 
punch, and all the addition and formation of differences and 
summations are done on the tabulator. In all this work with 
the tabulator we assign the first counter to the argument and 
to the identification of the particular step in the calculation, and 
counters (2), (3), (4) to the corresponding numbers pertaining 
respectively to x, y, z. The work of the operator consists of 
passing cards from one machine to the other, the selection of 
a card from the files of tabular functions, and the inspection 
of the highest order differences printed. The calculations 
are always repeated in a given sequence which is quickly 
learned. All switching of the machines to change from one 
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operation to another is performed automatically by means of 


_ the calculation control switch. Thus the machines are always 


ready to perform the next computation when the last has 
been completed. 

The extrapolation of X, X”, --- for equation (2) requires 
only a few seconds, and the same is true of equation (3). 
Take for example the extrapolation of X one step in advance. 
It is only necessary to add together all the differences down 
the diagonal from the right to left. This gives the extra- 
polated quantity with an error of the lowest difference neg- 
lected. In practice where fifth and sixth differences are em- 
ployed the highest should always be less than one hundred. 
The factor 1/12 as well as the small factor contained in the 
function X should be sufficient to give 1/12 X correctly. 

The calculation control switch contains a row of electric 
contacts each of which is operated by a rotating cam. The 
cam is a circular fiber disk which is notched at various points 
around the circumference. A series of about twenty of these 
disks are attached to a common shaft to form a sort of player 
piano roll. When this roll is rotated from one position to 
the next the various contacts open and close according to the 
notches in the disks. The circuits from the contacts are used 
to operate the various control switches on the tabulator and 
multiplier, and also a number of multicontact relays which 
effectively change the wiring of the plugboards. Each step 
in the integration consists of a certain number of distinct 
machine operations which always come in the same order. 
Hence in order to have the machines ready for each operation 
it is only necessary to rotate the roll from one position to the 
next, one complete revolution corresponding to a complete 
step in the integration. One roll serves for all equations of a 
given form, and a new one could be prepared in a few hours. 

As an example of the method consider the integration of 
the equations of planetary motion, which has been performed 
for many different planets. A detailed discussion of this 
problem will be given in chapter XII, but we shall describe 
it here in general terms. 
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For planetary motion equation (1) becomes 


P= 


x; = a function of ¢ alone, 

p2? = (x; — x)? + Oi — 9)? + (2: — 2)%, 
m; = small constants, 
P=xer®tyt+s. 


Since m,; are small we use approximate values of x, y, s and 
compute P in advance. P, for our purposes, is thus a func- 
tion of ¢ alone which is computed and punched on a card for 
each value of the argument. The equation then becomes 


kx 
—— at Pit). 


A table with argument 7° and giving Z is also made up in 
advance. This works for all planets. During the process of 


integration, r° is computed and sis found by interpolation in 
the table. 

For this integration ten significant figures are used in the 
calculations, about 60 points are used for each revolution of 
the planet, and fifth and sixth differences employed. The 
integration proceeds at the rate of about 12 points per hour. 

It might be added that even if the tabulator and summary 
punch only were available, it might in some cases be efficient 
to evaluate X with an ordinary electrical computing machine, 
and use the tabulator for all the differences, summations, and 
estimations. 

The Punched Card Method may be used in connection 
with other methods of solving differential equations. For 


example, chapter XI illustrates solution by the substitution 
of harmonic series, 


PART II]. ASTRONOMICAL APPLICATIONS 


CHAPTER IX 


THE CALCULATIONS INVOLVED IN THE 
CONSTRUCTION OF A CATALOGUE OF 
PHOTOGRAPHIC STAR POSITIONS 


1. INTRODUCTION _ 

The most important catalogue of this type is that being 
constructed by Dr. Frank Schlesinger and Dr. Ida Barney at 
the Yale Observatory, and known as the Yale Zone Catalogue. 
It is described here because of the part played by the punched 
card method in performing the laborious calculations which 
are inherent in the problem. The punched card method was 
not adopted until the program had been under way for several 
years, and for this reason the methods previously employed 
were changed only where the advantage of the new method 
was very great. Had the punched card method been avail- 
able at the outset it would no doubt have been used to fuller 
advantage. As it is, a saving of fifty per cent has been intro- 
duced in the computations which are done on the cards, and 
as these constitute about half the complete computations, a 
saving of twenty-five per cent has been introduced in the total 
cost of the calculations. The calculations now performed on 
the cards illustrate many of the methods of part II of this book. 

The Yale Zone program consists of the reobservation of 
the relative positions of all the stars in the so-called A. G. 
Catalogue. This catalogue was completed about forty years 
ago, so the comparison of the old and new positions gives very 
accurate values for the relative motions. The positions in 
the A. G. Catalogue were obtained by visual measurements 
with a meridian circle and constitute the combined work of 
several observatories for many years. The photographic tech- 
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nique which is being used in the new program is so efficient 
that the calculations now become an important part of the 
total cost of the project. 

In the photographic method a series of photographs of the 
sky are taken, and the positions of the stars on the photo- 
graphic plates are measured with a measuring engine. From 
these measures of rectangular coordinates (X, Y) on the plate 
must be computed the spherical coordinates on the sky, 
These spherical coordinates are known as right ascension and 
declination (a, 5). The computation must be made for each 
star as it appears on two overlapping plates, and with seven 
significant figures. In order to facilitate the identification of 
the stars in the measuring engine the process is reversed, i.e., 
the spherical coordinates given in the A. G. Catalogue are 
converted into the rectangular coordinates as they appear on 
the plates. Then, since the relative motions of the stars in 
fifty years are minute, the small differences between the ob- 
served and computed X and Y may easily be converted into 
the small differences in a and 6. The application of these 
small corrections to the old @ and 6, give the new a@ and 6, and 
the long conversion does not have to be made again. The 
most costly calculation in the process is thus the conversion 
from spherical to rectangular coordinates. The well-known 
formulae are: 

tan Aw 
(1) + © eer acieen ft cos 8’ 
tan 6’ sec Aa cos 6 — sin 6 


Y= hae co daaa eT cost 


X, Y are the rectangular coordinates of a star on the plate, 

6’ is the declination of the star, 

§ is the declination of the center of the plate, 

Aa is the difference in right ascension between the star and 
the center of the plate, and 

K is a constant for the entire program. 


The average number of stars on a plate is about 1,000 and the 
plates overlap so that each star appears on two plates. 
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Special tables have been printed for effecting this trans- 
formation but when the plates are large, and accurate posi- 
tions are required, the use of the tables becomes more expen- 
sive than the direct calculation. 

The second computation which is now done on the cards 
is the computation of the precession for each star. Right 
ascension and declination are measured with respect to a 
moving coordinate system, so that the changes in these 
coordinates consist not only of the small relative motions of 
the stars themselves but also of the motion of the axes which 
is much larger. In order to find the so-called “proper mo- 
tions’’ it is therefore necessary to allow for the “ precession.” 
For the purposes of the present catalogue the precession in a 
and 6 for a given star may be computed from the following 
formulae: 

(2) a = ao + pat + Sal”, 

6 = 69 + pat +- Sif’, 

ao, 69 = right ascension and declination to be corrected, 

a, 6 = corrected right ascension and declination, 

¢t = number of years for which precession is to be applied. 
In this case ¢ is the same for a very large number of 
stars. Pa, Ps, Sa, Ss are functions of ao, 69. They have 
been computed for each star and are given in the A. 
G. Catalogue. 


The computation is simple but must be made a great many 
times. It is complicated by the fact that astronomers still 
use the sexagesimal system for star positions. a, a are given 
in hours, minutes and seconds while 6, 59 are given in degrees, 
minutes and seconds. 

In order to make the results of such an observing program 
immediately available the observations are made in zones of 
declination. Inthe A. G. Catalogue, for instance, the various 
zones were constructed at different observatories. The new 
photographic catalogue is all being prepared at Yale, but the 
various zones are published as soon as they are completed. 
The following shows the work done in this program and the 
part played by the punched cards: 


Plates 
Zone Taken 


+55 to +60 1916 1926 1930 
+50 to +55 1916 1916, 22 1925 
+25 to +30 1928 1928 1933 
+20 to +25 1928 1928 1934 
—2to+ i 1914 1914 1926 


Measured | Published | Punched Card Computation 


+ 9to +20 1940 1940 Coordinates, Precession 
— 2to+ 9 1936 1939 Coordinates, Precession 
— 2to —10 1933 1938 Coordinates, Precession 
—10to -—i4 1933 1934 Coordinates 
—i4 to —20 1933 1934 Coordinates 
—20 to —30 1933 1934 Coordinates, Precession 


2. THE PUNCHING OF THE CARDS 


For each star in the A. G. Catalogue a card is punched 
containing all the pertinent information: the star number; 
the right ascension in hours, minutes, seconds and hundredths; 
the declination in degrees, minutes, seconds and tenths; the 
precession and secular variation in a and 6; the mean epoch, 
and the magnitude. The last two items are for a purpose 
other than the computations here described. In addition to 
the customary verification of the cards, totals have been 
formed for some of the columns of the catalogues and these 
totals used for checking. 


3. THE COMPUTATION OF THE PRECESSION 
The precession for an interval of fifty years is to be applied 
according to formulae (2). The necessary calculations for 
each star are: 
a= a+ 50. + 12.5s,, 
= 60 ss 50 Ds se 12.55;, 
where ay, 50, Pa, Ps, Sa, $3 are given on each card and 
292 < pa < 338, 
ps < 2071, 
Sa < 090250, 
S; < 07600. 
The results are required to 0°001 in a and O01 in 6. 
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The necessary multiplications, additions, and subtractions 
could be made by successive operations on the multiplier, 
but for such a large number of cards it is much cheaper to 
perform part of the work on the sorter and gang punch. By 
this method about 98 per cent of the reductions.in a, and 80 
per cent of the reductions in 6 require only one run of the 
multiplier. The secular terms are computed with the aid 
of the sorter, and punched with the gang punch as true figures 
or as complements. The work on the multiplier is always of 
the form 

A+t+BXC+D, 


where A is the uncorrected coordinate, B X C is the precession 
and D is the secular term which has been applied with the 
gang punch. D also includes such additional digits as are 
necessary to convert excess minutes and seconds in the 
s2xagesimal system. These factors must be determined be- 
fore the operation on the multiplier. For most of the stars 
this is easily done with the aid of the sorter, and for the 
comparatively small remainder, a trial run on the multiplier 
is necessary. 

In the reduction of the right ascension the precession is 
always positive. Hence if positive secular terms are punched 
as true figures and negative ones as complements the operation 
on the multiplier will always be the standard one 


A+(BX CO 4D. 


The secular variations range from — 0.0250 to + 0.0250, 
which means that the cards must be divided into 500 groups 
of about 60 cards each. This requires three sortings. The 
gang punching is most easily done by taking every eighth 
group in order that 12.5 sq will change by one unit. 

Before starting the multiplier the cards must be divided 
into groups according to whether 60, 120, or 180 seconds are 
to be converted into minutes, and 0 or 60 minutes are to be 
converted into hours. Let us consider first the question of 
the seconds. The secular term is never greater than about 
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083, and hence can affect only those stars which are near the 
critical values. In making the separation to the nearest 110, 
or for about 98 per cent of the stars, we can ignore this 
quantity. In order to save time at the sorter we make this 
separation in two steps, first making the separation to the near- 

_est 10: and then to 1s. In each case we first divide the cards 
according to the number of seconds in a» and then for each 
group find the critical value of pe. 


Convert 0° 


Fic, 14, Computation of Precession in Right Ascension 


The sorting and gang punching is simplified by the use of 
the diagram shown in figure 14. We could plot each star 
on this diagram according to the seconds of a and the value 
“ ba. Those lying below and to the right of the shaded 
diagonal must convert 60 into 1", while those lying to the 
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left and above must not. For those within the shaded area 
it cannot be decided without considering the value of 
12.5 sa. If they are separated according to the sign of s. 
half of these can be eliminated. There are also two tiny 
regions of uncertainty in the corners due to s,. It is obvious 
that for the majority of the stars the decision can be made by 
sorting only in the tens of seconds in a» and the tenths of 
seconds in fz. A large portion of the remainder requires 
sorting on the seconds of a» and the hundredths of p,. The 
process is continued to the point where it is cheaper to make 
a trial run on the multiplier than to continue the sorting 
process with small groups. 

The trial run on the multiplier is necessary only for a small 


60 + &e = 
P§ +)ors Ps — 
$6 — $65_-* 

8 "0 1-2 


Convert 60° 


Fic, 15, Computation of Precession in Declination 


We bya nee 
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percentage of the cards. In order to avoid punching these 
few cards and thereby making them different from the others, 
the operator simply watches the numbers as they appear in 
the summary products counter and separates the cards as they 
come out of the machine. 

As soon as the number of seconds to be converted is de- 
termined we can sort the cards on the minutes of a and 
determine which ones are to have a carry-over into the hours. 

Since we now know in advance where all the sexagesimal 
conversions are to be made, we can incorporate the necessary 
quantities into D of the machine operation A + (BX C)+D. 
The punched result is thus the right ascension corrected for 
p. and s, and expressed in the customary system of hours, 
minutes, seconds, and decimals. 

For the precession in declination the problem is compli- 
cated by the fact that 4) and ~s; may be positive or negative 
and by the fact that s; is larger than s.. The cards are there- 
fore sorted first according to the signs, and separate diagrams 
are used for the different cases. The diagram in figure 15 is 
the one used for the case 69 +, ~; +, 5; —, or for the case 
5) —, Ps —, Ss +. In this diagram there are the added di- 
agonal lines indicating values of 12.5s; = 2”, 4’, 6”, 8’. In 
this case also the test run on the multiplier is used for the 
doubtful cards. 


4. THE COMPUTATION OF THE RECTANGULAR COORDINATES 


This requires the evaluation of formulae (1) for each star 
and for two overlapping plates, with an accuracy of seven 
significant figures. This calculation illustrates the extensive 
use of special tables in the form of files of punched cards 
(see ch. IV). The special tables used are: (1) a table of 
tangents with the tabular interval 10”. This table was made 
from the values in Brandenburg’s table. (2) A special table 
of tangents and secants for angles less than about 6° and with 
a tabular interval of 0!1. This table was constructed for the 
purpose with the aid of Brandenburg’s table of cosines. 
(3) A table of reciprocals multiplied by K, for all values of 
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the denominator in formula (1). This table was constructed 
for the purpose. 

The features of the machine work consist mostly of avoid- 
ing multiplication where possible, and of going from one card 
to another ‘at the proper time. 

The first card which was punched from the A. G. Catalogue 
is used also for the computation of the precession. There is 
still some space left on card (1) and this is used for tan 6’ and 
the tabular difference. There is not room for the interpola- 
tion on this card, but having found the tangent and difference 
it will not be necessary to carry all of 4’ itself to the next card. 
The application of the tangent and difference is of course 
accomplished automatically by the method of chapter IV, 
section 3. The cards are tabulated for a check. 

A second card is now punched with the high-speed repro- 
ducer, transferring from the first card the star number and 
the minutes and seconds of the right ascension. At the same 
time the gang punch applies the complement of the right 
ascension of the plate center, with an additional factor of 
3600 for those stars which require a carry over of one hour. 
It should be remarked that the right ascensions of the plate 
centers have been rounded to the nearest whole second so the 
decimals will not have to be subtracted. 

One run on the multiplier of card (2) now gives as a cross- 
footing operation Ae; and Aas in seconds of time. Checking 
is done on the tabulator. 

The values of the tangent and of the secant of Aa and of 
Aa, together with the differences of the secant are now applied 
and checked in the usual manner. The differences of the 
tangent are not necessary, as they may be inferred from those 
of the secant. 

The card files of tan Aa, sec Aa have been made with a 
small interval of tabulation to facilitate interpolation. This 
was done because the range of argument is small and there are 
four interpolations for each star. The tabular interval is 0:1, 
so in the seventh place the difference in the tangent is about 
seventy-five and varies very little, while the difference in the 
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The numerator of the Y equation is next to be computed. 
The value of sin 6 is gang punched on the card and the cosine 
comes from a rate card. The numerator is more troublesome 
_ than the denominator for zones near the equator since cos 6 
is larger than sin 6. We therefore have to differentiate be- 


: 1 
secant is never more than eight but varies considerably. The | 
| 
7 
tween cases where tan 6 = tan 6’ sec Aw. Checking is done i 
{ 
i 


interpolation of both is done by one crossfooting operation. 
A set of master cards is sorted into the detail cards according 
to the value of the difference of the secant and according to 
the hundredth and thousandth of a second of Aw. In addition 
to the quantities used for sorting, these cards contain the 
number to be applied to the tangent and secant for the 
interpolation. The operation on the multiplier is done as rate 
card crossfooting. Thus with two sortings, one crossfooting 
operation, and one check run on the tabulator we are able to 


with the tabulator. 

From the special table we find K times the reciprocal of 
the denominator and interpolate the result. Checking is done a 
by progressive digiting. 


perform and check two interpolations. The tangents were | We now reproduce from card (2) on to cards (4) and (3) 
made by subtabulation from Brandenburg; and to avoid , ae the values of tan Aa, namo and sai ready = aa : 
accumulation of error, were taken to eight places, and are Bi and x by straight multiplication. ese are eight-dig + 
interpolated with eight. In order to save space on the card a) multiplications and require eight runs for each star, including \ 


the first three digits of the secant which are constant (1.00) 
are not punched. The above is done for Aq; and then for Aa». , 
The third card is now required. From card (2) are trans- is 

ferred the star number, sec Aa, and sec Aas; and with the | 
gang punch is added the plate number and the sequence of 

digits 1 23456789. From the first card is transferred Nils 
tan 8’ and its difference. Bia ties 

The tangent of 4’ is interpolated on the multiplier as an 
operation of the form A + B X C. Before this interpolation, 
the cards are sorted according to the difference so as to form 
constant multiplier groups which can be checked on the Ep ae 
tabulator. - Rae 

The product tan 6’-sec Aa is now punched and checked on : 
the multiplier. The first digit, “1,” of sec Aa which was not 7 
punched on the cards at the time of the interpolation is now Fi 
obtained from the sequence of digits mentioned above. 

The denominator of formula (1) is computed on the multi- | 
plier as an operation A + B x C, and checked on the tabu- iy hi 
lator. Sin 6 is taken from a rate card and cos 6 is made up ; 
from the sequence of digits. This is done for Aa, and then By 
for Aas. A 

Card (4) is now prepared. From card (3) comes the star 
number, the denominator, and the quantity tan 6’ sec Aaj. 
Card (5) contains the same quantities for Aa». 


the checking. | 
Finally the results are listed with the tabulator. 


ef 
| 


CHAPTER X 


STELLAR PHOTOMETRY 


1. INTRODUCTION 


The Rutherfurd Observatory of Columbia University is 
engaged in the determination of the relative brightness 
of stars by means of photography. The present program 
which should be completed in a few years includes about 
150,000 stars. While this work was started before the appli- 
cation of the punched card, it is very doubtful if the program 
in its extended form could have been undertaken had it not 
been possible to perform the necessary calculations by the 
punched card method. The process is now so automatic 
that the only personal attention given the individual star is 
the actual reading of the photometer. 

In the photographic method of determining the bright- 
ness or magnitude of stars two exposures are taken on each 
plate: one of a standard region containing a selected sequence 
of “comparison stars” and the other containing the stars 
whose magnitudes are to be determined, the so-called “pro- 
gram stars.” The plate is then placed in a photometer which 
measures the amount of light obscured by each of the star 
images, the amount of obscuration being indicated by the 
deflection of a galyanometer. The readings of the galvanom- 
eter for the comparison stars of known magnitudes give a 
calibration curve for converting the galvanometer readings 
of the program stars into magnitudes. The photographs are 
taken with a grating over the lens of the camera so that each 
star gives a central image and two diffraction images, and 
each star is photographed on two overlapping plates. For 
each star on each plate are recorded the reading of the gal- 
vanometer on the general background in the neighborhood, 
the reading on the central image, and the mean of the readings 
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of the two diffraction images. For the brightest stars the 


diffraction images only are used, and for the faintest ones only 


the central image. 

The readings on the star images must be corrected for the 
transparency of the background. The reduced reading for 
the star image (central or diffraction) is given by subtracting 
the reading on the star from the reading on the background 
and multiplying the difference by the ratio of the standard 
background reading to the observed. In the program here 
described the galvanometer readings range from 0 to 200, 
and 200 is the standard background reading. Hence if S is 
the reading on the star and G the reading on the background 
near the star the corrected reading is (G — S) a 

The corrected readings for the comparison stars may be 
plotted against the corresponding magnitudes to give a cali- 
bration curve for the plate. From this curve may then be 
read the magnitudes corresponding to the corrected reading 
for each program star. However, a study of a large number 
of these calibration curves has shown that within the experi- 
mental error they may all be reduced to a standard curve by 
a linear reduction. Thus if m is the magnitude read from the 
correct curve for a given plate and m’ is the magnitude read 
from the standard curve, then 

m =a +t bm’, 

where a and 0 are constants for the plate in question. This 
fact makes it convenient to substitute a numerical process 
for the plotting and reading of the graphs. The standard 
graph is made into a table from which m’ is read. The 
forty or so comparison stars are used to determine the con- 
stants a, b and their probable errors. The magnitudes of the 
program stars are taken from the table and reduced with the 
a, b obtained from the comparison stars. 

In addition to the above reduction of the observations, 
there is the problem of identifying the stars on the plate for 
measurement. This identification is made by means of rec- 
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tangular coordinates which must be com 

! puted from the know 

spherical coordinates. on 
When the magnitudes have been determined for a large 

number of stars they are examined for systematic errors de- 

pending upon the magnitude, the position on the plate etc 


2. Tue Use or THE PUNCHED CARDs 


The first step is naturally the punching of the cards fro 
catalogues such as the B.D. The data for each star 
punched on a card and verified. Each card contains he 
number of the zone in question, the number of the star in th 
zone, the right ascension and declination, and the visu ‘ 
magnitude. With the gang punch is now added the'n be 
of the plate on which the star will appear. ha 

The computation of the rectangular coordinates on the 
sa is the same problem as that of chapter IX except that 
els seal purpose very little accuracy is required. 3a 
pene seen the series expansions in terms of Aa, Aé are 
eat ante = yond the second order are neglected. To 
chosen that pared vail b. oS pian 
: : ey all have the same declina- 
ses and the difference in right ascension between ie, i 

nstant number of whole minutes. Further, the scales 


on the photometer in x 
and y were ruled in mi 
The formulae required are: se andes ot arc. 


A= kA — k2Aa’ASs’, 
Y = Ad’ + k;(Aa*)?, 
ki, ke, ks are constants for the zone, 


Aa* and Ad’ a i i 
re the differences in a and 5 between the star and 


th 
; c center of the plate, the first expressed in seconds of 
e and the latter in minutes of arc. 


The : 
a star eee ae of Xi, V1 and Xs, Y2, the coordinates of 
amount of sorti lerent plates, is accomplished by a small 
ing and gang punching, the pulling of one card 
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from a file, and one summary punch operation. The first 
operation with the sorter and gang punch gives 46’. The 
second gives the minutes of Aa. The third gives the small 
term k,Aa'Ad’ for one of the plate centers, and then for the 
other. A card is now drawn from the file with argument 
Aa™* and giving k\Aa* and k,(Aa")’. This card is still taken 
from the file by hand because there are very few of them and 
because the question of duplicates is troublesome when the 
sorter is used instead of the collator. The tabulation of the 
original card and that pulled from the file gives at one run 
the coordinates of the star as it will appear on two plates. 
The summary punch gives a new card containing the zone, the 
plate number, the star number, the visual magnitude, X;, ¥1, 
Xo, Y2. 

The cards just obtained are to be used at the photometer. 
To facilitate the reading of the cards they are passed through 
the interpreter which prints at the top of the card the numbers 
which are punched on the card. Thus, in type, at the top of 
the cards appear the number, the visual magnitude, and the 
four coordinates. 

The observer uses the printed numbers to locate the star, 
and then places the card ina card punch to record the readings 
of the galvanometer directly on the cards. By the time the 
galvanometer has come to rest the card is in position for 
punching, and the scale reading for background, central image, 
and mean of the diffraction images are recorded. Each card 
is used at the photometer in connection with two plates. 

For the machine application of the background correction 


we rearrange the formula: 
200 


Corrected reading = 200 — | -S. 


ames : . 1!) 1 ee ; 
The instrument is adjusted so that — varies little. The 


200 . d ie 
cards are sorted on G, and G's used as a rate card on the 


Since G differs very little from 


Re 200... 
multiplier to compute —— -S. 
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200 we can form this quantity for both the central and the 
diffraction images with one run of the machine. This is done 
as an operation A + B X C where A is the reading on the 
central image and the reading on the diffraction image, B is 
200 : 

Tr 1, and C is the two largest digits of each of the two 
parts of A. As an example take the case where the reading 
on the background is 195, the reading on the central image is 
136, and the mean of the diffraction images is 183. Then 


A = 13600183 
B .026 
C = 13000018 


The answer is 139.38, 187.68 which are rounded to 139.4 and 
187.7. These quantities may be in error by as much as 0.26 
because of the omission of the end digits of C, and by greater 
amounts for larger values of 200 — G. 

The table giving the provisional magnitude m’ is made in 
the form of a file of punched cards which may be used in the 


: 7 
usual manner. To avoid the subtraction of -~ -S from 200, 


the 200 is absorbed in the argument of the table. The opera- 
tion of “looking up” the number in the table requires about 
five hours for 20,000 entries. 

The values of m’ for the comparison stars must now be 
used in conjunction with the corresponding values of m for 
the determination of the plate constants a and b by the 


method of least squares. There are m (about forty) obser- 
vation equations of the form : 


a+m'b—m = 0. 
The normal equations are 


n-a +([m’]-b —[m] = 0, 
[mn'}-a + [m'?]-b — [mm’] = 0, 
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and the solution is given by 
[om }[m'* | — [im ]mm'] 

n{m?|—[mP ’ 

wr de n{mm’ | — [m][m'] 
n{m’ * —[m'P ° 

The formation of the coefficients of the normal equations is 
done very simply on the machines. The card for each com- 
parison star contains m and m’, so mm’ and m’ can be 
formed with two runs of the multiplier. (Two runs are re- 
quired because of the wide variation in a and b.) The tabu- 
lator then gives with one run of the cards n, [m’], [m], [m'’?], 
[mm }. If relatively few plates are being reduced at one 
time the normal equations can be solved by hand from the 
quantities given on the tabulator record, and summary cards 
are not required. If say five hundred plates are reduced 
at once it might be desirable to eliminate this handwork. 
The summary cards are run through the multiplier to give 
[m ]Lm’*], (im! ][mm'], n[m’], (m'}, n{mm'], [m][m'). 
The crossfooting device then gives the denominator and the 
two numerators. A small table of reciprocals is required to 
avoid the division. One run on the multiplier for the inter- 
polation of the reciprocal and two more for the multiplication 
are then required to give @ and b. 

Having found a and } for each plate, we can evaluate m 
for the program stars by a run on the multiplier of the form 


A+BXC, 


where A =a, B=), C= m'. This quantity should be 
the comparison stars in order to 


formed at the same time for orde! 
determine the weight of the solution fora and 6. This weight 


is determined from the sum of the squares of the residuals 
where the residual v = m — 2 — bm’. 

The values of m may now be used to form the mean of the 
observations for each star and to form the difference for each 
star of the values obtained from the two plates. These are 
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crossfooting operations except for the multiplication by 4 
in forming the mean. 1 is here 2, 3, or 4 depending =e 


whether or not central and diffraction images were measured 
on both plates. The values of the difference for the he 
plates are used to discuss the accuracy of the results The 
values of the mean of this quantity for a given plate mi b. 
tracted from the individual values, and the outstanding ng 
crepancies are reobserved. These outstanding differences 
sometimes show that the wrong star was observed on one f 
the plates or that the star is variable. Finally the differen : 
are examined for systematic effects depending upon the ia 
nitude, the X and Y coordinates on the plate, and the ie 
(difference between photographic and visual magnitude) All 
aa ca are on the card and the necessary computa- 
i phe one most effectively on the tabulator and counting 
Jie de salerangall magnitudes are listed with the tabu- 


CHAPTER XI 


NUMERICAL LUNAR THEORY 


1. INTRODUCTION 


In this chapter the solution of a set of differential equa- 
tions is determined in the form of harmonic series with 
numerical coefficients and with arguments which are linear 
functions of the independent variable. An approximate solu- 
tion is substituted in the equations, and from the residuals 
corrections are obtained for the coefficients of the assumed 
solution. The project here described is the most extensive 
which has been done, and it was feasible because of the avail- 
ability of the punched card method. 

It is well known that the solution of the main problem of 
the lunar theory may be expressed with great accuracy as 
harmonic series giving the coordinates in terms of the time. 
The argument of each term is a linear function of the time 
with four fundamental frequencies. The coefficients in these 
series and two of the frequencies, the motions of the perigee 
and node, are determined from the differential equations of 
motion. In the well-known theories of Delaunay and Brown 
these coefficients and motions were determined by algebraic 
rs of the arbitrary constants, numerical 
at theend. These 
d no use is 


expansions in powe 
values of the constants being inserted only 
developments are built from the ground up, an 
made of previous developments except as a general guide. 
Hansen’s theory, though numerical, followed the general plan 
of the algebraic theories and was independent of other solu- 
tions. 

The basic principle involved in the present development 
is essentially different. It consists of substituting in the 
differential equations of motion the numerical values of the 
best available previous solution. The residuals thus obtained 
are then used with the equations ol 
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variation to determine 
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small corrections to each coefficient and to the motions of the 
perigee and node. This method gives at the same time a 
precise verification of the previous work, and increased accu- 
racy where necessary. It utilizes the previous theory, but the 
final results are independent of it. The coefficients involved 
are ordered solely according to numerical value and not 
according to the algebraic order of the parameters involved. 
The only assumption is that the square of the errors may be 
neglected. 

The numerical method by itself would leave much to be 
desired, but when used in connection with a good literal 
theory it gives results which could be obtained by the literal 
method alone only at the cost of tremendous effort. Thus 
in the present case the general accuracy of the theory of 
E. W. Brown was ample for observational purposes, but a 
few coefficients, particularly those giving the motions of the 
perigee and node, had not been obtained with sufficient accu- 
racy. To obtain the additional accuracy by an extension of 
the original method would require years of highly specialized 
work. Then too there was the possibility of error in the 
original work which had required twenty years to complete. 
By the numerical method it has been possible in about two 
years to verify completely Professor Brown's results and to 
obtain the required additional accuracy. 

The general method here employed was suggested by Sir 
George Airy about half a century ago. He spent fourteen 
years and considerable sums of money on the project but 
was unable to complete it successfully.* The numerical work 
involved is so extensive that with the methods then available 
numerical errors were almost unavoidable. In addition, as 
pointed out by Radau,j he had made one or two mistakes in 
theory. 

2. THE PROBLEM 


The problem investigated here, known as the main prob- 
lem of the lunar theory, is the determination of the motion of 


* Monthly Notices of the Royal Astronomical Society, 1888, vol. 49, p. 2. 
tT Bulletin Astronomique, 1887, vol. 4, p. 274 
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the moon under the gravitational action of the earth and the 
sun, with the assumptions that the three bodies are spherical 
and that the sun moves in a fixed ellipse about the center of 
gravity of the earth-moon system. 

The position of the moon is referred to rectangular axes 
through the earth’s center with x and y in the plane of the 
sun’s orbit, and rotating about the z axis with the mean angular 
velocity of the moon. ; : 

The equations of motion may be written in the homogene- 


ous form. 
(1) «Dy — yDx + (x? + y?) 


(2) 4D*%(r) + 40? + 9”) - 4[ (Dx)? + (Dy)? + (Dz)*] 
an ao aa) 108], 
= 2 +[-Z]+ p-| 2 (=; Yax J nat Sita 
7 oQ x ae 
(3) xD*s + | - De + x + 2Dy +3 — 32 
where x = 1 + x, and *, y, sare the deviations from circular 
motion in the plane of the orbit. 
Q = Lak jnwjn, 


i rers of as e St< and 
where k;, contains successive powers of a small constant a 


each w;, is the product of 
of order j by a harmonic se 
known numerical coefficients. 

The adopted solution is of the form 


a 4 


a homogeneous function of x, ¥, 2 
ries with argument ¢ and with 


- = La; cos ai, 
+ gee Db; sin aj, 
‘s = Ye; sin a, 
a; = (jp +ka t+ lr + ms)t 
j,k, l,m = 9, + 1+#2::- 
Ds q,.7,5 = given numbers. 


ed by multiplying the coefficient 
f tin a;, and changing the 
For each term, the 
To substitute 


The derivative series are form 
of each term by the coefficient of | 
signs where the harmonics are cosines. 
numerical value of D is the coefficient of /1n ai. 
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this solution in the equations (1), (2), (3) it is necessary to form 
all the product series of the second order for the parts in 
dependent of 2, and for the terms of second and higher rdler 
involved in. Since the individual parts of 2 contain ie 
factors, these latter are not required to the same Eecnrace wx 
the former. By far the largest part of the work consists of 
forming the product series for the second order terms. The 
coefficients in the assumed solution were taken with ey 
significant figures and there were about 500 terms for eact 
coordinate. aa 
For comparison with observation the rectangular coordi- 
nates must be converted into polars. The equations f 
this are: ie 


tan (V — nt) 


sin } 


aad again the transformation may be checked by substituting 
a numerical coefficients for the series representing x 
, ¢, 7. We find from the above equations 
(x* + y*)Dv = xDy — yDx, 
= ex + 1/6(sr)* + 3/40(2r)5, 


rx? = 1. 


> ¥. Ss 


3. THE : y z 
1E METHOD WITH THE PUNCHED CARDS 


The me >scri i f 
ated pi seit in chapter VII for large series was 
yed here. e multiplications were : ; 5 
ations were done as trig 
metric series i > i a tat 
eae iy = ep cases and as exponential in others, so that 
ote fi = Oo i at chapter occurred. Conjugate imaginaries 
=x-+ vy and s =x — ) 
pee, 8 5 = x — wy were employed to facilitate the 
ieee - of 2 and its derivatives. There were about eight 
dred terms in u, and i i | 
oy: in forming 1? and al 
x TEN dara and us about twenty 
psa d as products greater than 1 X 10-" had to be 
siderec 2g Z ' iff: 
or the sum and for the difference of the arguments 


CHAPTER XII 


THE COMPUTATION OF PLANETARY 
PERTURBATIONS 


1. THE PROBLEM 


The computation of planetary perturbations furnishes 
excellent examples of the methods of chapters VI, VII, and 
VIII. The methods of chapters VI or VII or both would be 
used in the so-called method of general perturbations while 
the application of chapter VIII is known as the method of 


special perturbations. 
The differential equations of planetary motion are 
a , kx 
aig ta a 
dt* ? 
with similar equations for y and 2. 
oy 
i ee , 
pi 
2 = (x; — x)? + (yi — 9) t+ 
, = ratio of the mass of the distur 
the sun. 


bing planet to that of 


The integration of these equations is always based on the 
fact that m; is small, never greater than 1 X 107°. When 
m; = 0 the solution is an ellipse with the well-known law of 
areas. The elliptic solution may be substituted in the right 
hand side of the equation and a second approximation to the 
solution made which gives the so-called first order pertur- 
bations. The substitution of this 
order, etc. 

The solution of the complete equat 


the form of a table of numerical values of x, 
through a certain range, 
arguments are linear 


then gives the second 


ions may be given in 
y and z for equi- 
distant intervals of time or in the 


form of harmonic series with terms whose 
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functions of the time. The expression of the values of the 
coordinates as a table is the method of special perturbations, 
while that employing the harmonic series constitutes the 
method of general perturbations. The advantage of the 
method of general perturbations is that it gives explicit ex- 
pressions for the coordinates in terms of the time. They may 
be used for computing the position of the planet at widely 
separated times without computing all the intermediate 
positions. The determination of these expressions with high 
accuracy is however a very laborious process. The table of 
numerical values is useful only for the range of the table, but 
it has the advantage that for limited intervals of time the 
table may be constructed with great accuracy and com- 
paratively little work. 


2. THE COMPUTATION OF SPECIAL PERTURBATIONS * 


This problem furnishes an ideal example of the method of 
chapter VIII. There are nine planets with m; large enough 
to be considered and about 1,500 planets and about 500 comets 
to which the process could be applied without changing the 
machine setup. In other words the necessary initial cost of 
the application of the punched card method can be distributed 
over a great many Cases. 

The card files required are: (1) a table of quantities to 
the minus three halves power and (2) tables of the coordinates 
of the nine perturbing planets. The first of these tables was 
made from Comrie’s table and is accurate with linear inter- 
polation to eight places. The coordinates of the perturbing 
planets were prepared from Comrie’s table and other sources. 
The coordinates are tabulated at ten-day intervals with five 
decimals in x, y, z. In order to use the center of the sun as 
origin instead of the center of gravity of the solar system we 
have on the same card the attractions of the perturbing 
planets on the sun. 


Th e ee = Xi = X. 
rhe quantity Mim is called the perturbation in x, 
, 


et The Computation of Special Perturbations by the Punched Card Method. 
W. J. Eckert. Astronomical Journal No, 1034, 1935, 
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and there are similar expressions in y and z. These perturba- 
tions are computed before the main process of integration is 
started. Since m; are small quantities we can use approxi- 
mate values of x for this purpose. Obviously if elliptic 
values are used for x, 9, 2 the integration will give coordinates 
which include the first-order perturbations. if first-order 
perturbations are included in the approximate coordinates, 
the results will contain second-order perturbations. 

The computation of the perturbations isa straightforward 
process on the machines. It is most efficiently done as mass 
production. The more disturbed and disturbing planets ~~ 
the more dates included the more efficient it becomes. 
cards containing the approximate values of x, y Z are = 
together with the cards containing the oo, fi fe 
and the sorting is done according to the date. . : si a 
and summary punch then give a new card for each con: 
taining (x; — *), (Wi — y), (2: — 2)- Repetition of ; wha 
ess gives successively the corresponding ge cae ia 
3, ---, 9, or as many of them as are required. : ie ca nh ie 
carry in addition to the date and the ne 
ordinates, numbers to represent the values of ¢ for the disturb 
ing planet and of j for the disturbed planet. Sede eh es 

When all the differences ol a baie EscnIee 
they are passed successive y throug! - multipy 
|. ei x)?, (yi — 9)": (zs; — 2)? and pi. | a ss 
it is necessary to make sure that no accuracy 1s los 


ne disturbing and disturbed planets, 


close approaches of tl pgs 


since the force increases sa 
method for doing this is given in th 


as the distance decreases. 
e last reference. 


to the values of p?* are 


1 
The values of —; corresponding 
pi 


1er with the sorter and 
is done with the multi- 
on the multiplier. 


taken from the table in the usual mant 
high-speed reproducer. proapolation | 
plier. The multiplication by ™m: is ¢ one yense pes 
In addition there is a factor (104)-*? which n cis st 
because our table of quantities to the minus 


> fr 2 ) 20 whereas 
pow I S Stl a ar nts from Zt 
fer Was const ucted W ith argumec 


p2 is frequently outside this range. 
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Final multiplication by (x: — x), (vi — 9), (i — 2) gives 
the perturbations due to each planet. 

Since the computation of perturbations is complicated, 
they are computed for fewer dates when possible and the miss- 
ing ones found by interpolation. The subtabulation is usually 
to halves and quarters and is done automatically by the 
Lagrangian method described in chapter V section 3, 

Finally all the perturbations for a given disturbed planet 
for a given date are collected and added with the tabulator. 
The summary card is made in red to facilitate its identification 
during the process of integration. These perturbations may 

now be considered as a function of the time for the process 
of integration. 


The equations for the integration proper thus reduce to 


kis B5+PO 


the solution of which is given by 


x = "X + 1/12X — 1/240X” 4+ 1/1951X'Y + --. 
= "X + the “reduction,” 


with similar expressions for y and z. 
The following table shows the various cards used in the 
computation of the perturbations and in the integration. 


Information 
approximate coordinates, x, y, 2 
same with signs changed 
Xi — X, Vi — Y, Bi — 23 p?# 
computation of k*m,/ p*; 
multiplication of x; — x, y; — ¥, 2: 
sum of perturbations 


ith difference of the second summation of the total 
attraction 


extrapolated value of (42) 


integrated values of x, y, zand the values of the solar 
attraction computed from them. 


— 2, by k?m,/p’*; 


i v iat 
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Additional cards used in the interpolation of the perturbations 


are not shown. On each card appears the date in the form of 
the Julian day, the identification of the particular card in the 
general scheme, and the quantities required. Six columns 
are used for the date, and two for the identification. ards 
1, 2, and 3 are used in the computation of the perturbations 
while 4, 5, and 6 are used in the integration. 

The integration is best described by means of an example. 
Figure 16 is a portion of the tabulator record and shows all 
the totals printed during three steps of integration with 
seventh differences. Each line is printed at one time. The 
first group of ten columns comes from the first counter of the 
tabulator and contains the date and the identification of the 
operation ; the second, third, and fourth groups come from 
the second, third, and fourth counters and contain information 
pertaining tox, y; 2 respectively. In the first group of cr — 
digits 2 to 7 indicate the Julian day while columns L gee 8 
indicate the quantity contained in a given line. Thus from 
the seventh line (46) in the figure we find pete oh ae 
[27427980, X*Y = 24948, 0°" = + 1585, and 2 ie 
and find also from the line (40) that MI = 1.8462948 . i 

Summary cards are punched for each line pee ee 
identified as (49) and the second (40). The first sees é a 
are steps in the computation of the last (49) which contains 
the true seventh differences of the attraction. ae 

There is in addition a card with the same Sep ges 
as (60) except for a J a punched in columns -80. s 
card is used in computing the solar attraction. 

In the machine process we consider a step comp 
we have just finished computing the differences sae He 
tions. Thus on the tabulator record have eo? si ee 
lines (49), (48), °°» (40), and the corresponding em ad 
cards have been punched. The new step rope Sees 
difference and summation cards and proceeds gn 5 = 
tion of the quantities needed in the reduction. : stn 
(42) must be estimated one step advance * _ claire 
accuracy. This is easily done on the tabulator Dy > \ 
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the gravitational constant, and the result is ready for multi- 
plication by x, y, s. These quantities, if added to the per- 
turbations, would give the new values of (42) to be entered 
in the difference table. 

For simplicity of the machine operation we first find the 
highest-order difference and build the others including the 
two summations by progressive totals (see ch. V, sec. 4). 
The highest-order difference is obtained as follows. 


Vii = fi — [fo + Vo’ + vee Vor] 


where f is the attraction. In this case we add the solar attrac- 
tion and perturbation and subtract the differences for the last 
step. This leaves V,‘ in the counters of the tabulator, and 
by adding successively Vo, Vo'?, «++, Vom, Vo-? we obtain 
the values of all the required differences and summations of 
X, Y, Z. These quantities are printed on the record as 
(49), (48), ---, (40). One step of the integration is thus 
completed. 

The various operations in the above cycle are repeated 
for each step in the integration and always in the same order. 
All the necessary switching of the machine has been taken 
care of in the calculation control switch as will be described 
presently. The calculation control switch has twelve dis- 
tinct settings and these are run through in connection with 
the various steps outlined above. For six of the settings the 
switch moves automatically to the next, but for the others 
the operator must press a conveniently located button. 

We now describe the activities of the ope 
cycle. 

The calculation control switch is set on position 1. The 
operator takes cards (47) from summary punch, removes 
(40), (41) and adds a yellow master card and two blank cards. 
These are tabulated in the direct position. 
cards (52) appear while the switch automa 
positions 2, 3, 4, 


Operator takes a blank card, a salmon rate card, and 


rator during one 


Three summary 
tically moves to 


im ‘red perturbation card and three car 
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ceiites the following cards: (40) direct, cn rhe 
Bom multiplier reversed, the above cards from the : “i 
en aster, and three blanks direct. The summary pu : 
Bes four cards (60) while the switch moves automatically 
. 1 first three (60) cards in multiplier. While 
a peeretor hectic he lifts lever on summary punch = 
al ac and presses eject lever. He substitutes in = 
i of pete (60) card reversed for the blank. te t = 
. * ib in the summary products counter 0 in 
7 E attiol Dita operator pulls the proper card from 
q multiplier, 


e . . . 1 
“4 file of 3° 


r? are punched by hand. 


; ining digits of 
on (60), and the nave the switch to 8 and the 


3 a 3 iene eerie passed through the multiplier direct. 
3 ro 
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Move switch to 12, tabulate all (47) cards including (40). 
This completes one step in the integration. 

We have not indicated all the details of picking up and 
laying down the cards and putting them away. The efficiency 
of the process requires that these be done without lost motion. 

The following is the punching on the various cards: The 
ordinary detail cards in the hopper of the summary punch are 
prepunched with a “0” in column 71 and an “x” in column 76. 
The yellow master contains a “0” in column 10 and 5400012 
in columns 51-57 for forward integration and 5399992 for 
backward. The salmon rate card has 33333338 in columns 
41-48 and an “x” in column 46. The pink master contains 
“0” in column 10 and a “4” in column 51 and numbers in 
columns 52-57 which depend upon the backward or forward 
integration and upon the order of the highest difference. If 
7 is the order of the highest difference we have 40007 for 
forward and 3999; for backward. 

The summary cards (44) have four ten-digit fields: argu- 
ment and identification, functions of x, functions of y, func- 
tions of s. The summary cards (52), (60) have two ten-digit 
fields: argument and identification; function of x, or function 
of y, or function of s. The card which refers to x has “x” 
punches in columns 9-10, that which refers to y has “x” in 
column 9, and that which refers to z has no “x.” The ‘‘x’s”’ 
in these two columns operate two class selectors on the tabu- 
lator so that functions of x always go into counter (2), func- 
tions of y into (3), and functions of z into (4). In these 
operations with the summary punch there are multicontact 
relays which punch the successive cards from the proper 
counters and indicate the results with the proper combination 
OF “x's.” 

The tabulator and multiplier have had minor modifica- 
tions to permit remote control of their operation through the 
calculation control switch. These modifications become effec- 
tive only when the cables from the switch are plugged into 
the sockets on the machines, so that when the switch is un- 
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3 ine in the normal manner. The 
ac he machine operates in no 
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Switch Position 
45 6 7 8 9 101112/13 14 15 16 17 18 19 20 
Pace, ak xX xX 
x 


= SCMms Ane ene 
KXKKXKKXKXX XX 
XXXKXKXXXX 
MK x 


—_ 
~ 


wu re 
| tn de Gs 


-_ 
~ 


Not used 


@~ 


Cam No. 

1— multicontact relay no. 1 to select multiplicand 
2— multicontact relay no. 2 to select multiplier 
3— multicontact relay no. 3 to control punching of products 
4— multicontact relay no, 4 to control punching of products 
3s— multicontact relay no. 5 to control summary punching 
6— multicontact relay no. 6 to control summary punching 
i— multicontact relay no. 1A in tabulator 
s— multicontact relay no. 1B in tabulator 
9— multicontact relay no. 2A in tabulator 

10— counter clear on major control 

11— summary punch on minor control 

12— rate card multiplication 

13— A-BxcCc 

14— punch—check 

15— “x” punch in column 44, 
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3. THE Compuration or GENERAL PERTURBATIONS 
The method of general perturbations requires the expan- 
sion of the disturbing function 2 = and its derivatives 
+ 
as harmonic series with arguments which are linear functions 
of the independent variable. There is a wide variety in the 
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variables and of methods of expansion, and the 
“hice of the particular one to be used is dictated by the 
Sor f the problem. From the computational standpoint 
ver they are very similar, and the methods outlined here 
“san be readily applied to any case in point. | ; 
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-++, 11, and s = 1/2, 3/2, 5/2, and 7/2. It must be remem- 
bered that there is a different value of a for each disturbing 
planet. The required derivatives with respect to a are ob- 
tained from the same table, for 


aX (1 + a? — 2a cos S)* 


‘ x ees Dae (Ger — G,) cos iS, 
the factor 3 to be inserted for = 0. 

We thus see that even with the tables available the work 
of obtaining the necessary value of 8,“ and their derivatives 
is quite considerable. In order to perform this part of the 
calculation by the machines the tables have been punched on 
cards, each card containing a value of log G,“ together with 
the first and second differences. The selection of the proper 
cards from the file is easy, for there are only a few values of a 
required and all the cards for a given value of a should be 
kept together. The interpolation is done on the multiplier, 
and the factors (1 — a*)~*, a*, s are applied with the tabulator 
or multiplier. No provision has yet been made for finding 
the antilogarithms by machinery. This is a simple matter 
by hand with the ordinary tables, but it could easily be done 
automatically if sufficient use could be found for such a card 
catalogue. 

Once the series involving a have been found, the other 
expansions can be carried out by the methods of chapters VI 
and VII. Thus the elaborate and extensive calculations in- 
volved in the planetary theory may be done almost entirely 
without hand work. 
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aie seed to operate several sets of contacts (multi- 
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ADDITIONAL DESCRIPTION OF THE MACHINES 4 

AND PLUGBOARDS 


For those who may wish a more detailed knowledge of the 
machines we give here a few of the principles of their opera- 
tion. This is in no sense a complete description but is in- 
tended merely to remove some of the mystery. 

The regular tabulator service men are familiar with all the 
circuits of the machines and are glad to answer questions, but 
usually this is little help as the user of the machine is not in 
a position even to ask questions. It is the purpose of these q 
few paragraphs to give a basis for asking questions. 

The tabulator contains about 200 electromagnets. The 
passage of a current through one of these magnets causes an 
armature to be pulled toward the magnet, and the armature 
in turn is used to operate mechanical parts of the machine 
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Bureau has only four counters and five banks. In these cases 
we simply remember that some of the indicated holes have no 
meaning. 

We proceed now to describe the plugboard in terms of 
magnets and relays in order to make clear some features not 
easily grasped otherwise. In the tabulator most of the mag- 
nets and relays have one end of the primary circuit perma- 
nently grounded to the machine so that only the other end 
need be wired on the plugboard. In many cases two or more 
holes in the plugboard are wired in ‘“‘common”’ to the same 
circuit. This eliminates the need of double and triple plug- 
wires. 

At the top of the diagram are seven enclosures with ten 
pairs of holes in each. The two holes of a pair are common 
and are connected to the clutch magnets of the adding wheels 
of the counters and to the printing magnets. By wiring into 
these holes from the adding brushes it is possible to list and 
add at the same time. 

Directly under ‘No. 7 BANK” isa hole marked “cc.” This 
hole, wired to a position of the counter, gives the card count. 

A class selector is a multicontact relay with ten sets of 
three points. The diagram shows three class selectors. The 
primary circuits are each wired to a pair of common holes 
labeled “CLASS SEL. PICK-UP”’ and the thirty secondary con- 
tacts are wired to the groups of 30 holes labeled ‘‘cLAss 
SELECTION.” In using a class selector, say number one, the 
adding brushes from one field are connected to the ‘‘No-x”’ 
holes and those from the second are wired to the !‘x”’ holes. 
The “c” holes are then connected to the counter. The im- 
pulse to operate the primary comes from the ‘‘CONTROL 
BRUSHES.” 

The eighty holes labeled ‘‘ADDING BRUSHES” correspond 
to the eighty positions of the card. 

The next four rows of holes (“GROUP IND. & AUTO. CTRL.”’) 
are used for the control and group indication. For purposes 
of group indication there are multicontact relays with only 
two contacts for each circuit. When these points are closed 
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another hole in the plugboard which actuates the control 
mechanism. The failure of any one of the control relays to 
operate thus causes a “break in control.” By plugging into 
three of these positions we may secure the necessary impulses 
for major, intermediate, and minor control. 

These control relay points, if we skip the first, may be 
used for other purposes than control. When we require 
several relays in series which are operated by the comparison 
of holes in two cards, we simply wire in one of the holes and 
out another. (See part B of figure 21.) 

The next four rows on the plugboard are connected to the 
control brushes. 

Below these, on the right, are three pairs of common holes 
marked “Mayor,” “INTER,” and “Minor.” These are to be 
connected by plugwires to the holes above (‘‘To SHUNT PLUGS a) 
as described for the operation of the control. 

The holes “tcc,” “m(u),” “M(L)” are for special circuits 
used in listing cards with consecutive numbers to show that 
one of the sequence is missing. The hole “ucL”’ is for another 
special circuit. Our tabulator has neither of these features. 
The other holes in the section labeled ‘“x-pIsTRIBUTOR”’ are 
all associated with the operation of the counters. There are 
five double circuit relays for this purpose. These relaysare 
operated by ‘“‘x” punches on the cards. For each relay there 
are eight holes in the plugboard, the first two of which are 
common and go to the magnet of the relay. If one of these 
holes is wired to the control brush of a given column, an ‘x”’ 
in that column of the card will cause the relay to operate. 
The next three holes are wired to the three contacts of the 
first secondary circuit and the succeeding three to the other 
secondary circuit. When there is no “x” on the card the 
two circuits “‘c-No-x” are closed, but when there is an “x” 
the circuits “‘c-x" are closed. The five holes ‘‘PpLUG To Cc” 
supply electric impulses which are usually wired through the 
secondary circuits of the relays and used to operate the 
counters, 
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The last seven rows are used for the operation of the 

d list banks. The first column “App” is merely 
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change from the high tabulating speed to the lower listing 
bai! first group of switches shown on the chart is associ- 
ated with particular counters and list banks. The first seven 
switches are simple two-way switches. They have two de- 
scriptions in the chart: “TIST-IND.,”” and ‘“‘orr-ON.”’ In case 
there is a counter associated with the print bank the ‘“OFF-ON" 
applies, for in this case you can list or not list. If the bank 
is simply a list bank without a counter, provision is made to 
use this for additional group indication. In that case the 
switch causes it to list or indicate. 

The printing of totals from a given counter 1s controlled 
by a dial switch with five settings. The counter will print 
no totals; print totals on any break in control but not reset; 
print and reset on minor, intermediate, or major breaks in 
control. 

Negative numbers are added as complements in the 
counters. If the total to be printed is negative it may be 
printed as a complement, or as a true figure with a minus sign. 
The “BALANCE” switches are provided for this purpose. 

The “MIsc. SWITCHES” control the operation of the ma- 
chine asa whole. The “AUTO. START” causes the machine to 
start feeding cards after each reset. The ‘AUTO. RESET”’ 
gives a reset and total print operation on each break in control. 
The ‘AUTO. CONTROL” switches must be turned off when the 
control is not wired or there will be a break in control after 
each card. 

The ‘GROUP INDICATOR” switch operates the multicontact 
relay described in the explanation of group indication. This 
switch should be labeled “‘List-IND.” instead of ““OFF-ON.” 

The first summary punch switch, if “on,” prevents the 
tabulator from operating unless the summary punch is oper- 
ating properly. For instance the tabulator will not operate 
if the summary punch runs out of cards. 

The next three switches determine whether summary 
cards are to be punched on major, intermediate, or minor 
breaks in control. 
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The last switch is used mostly in bookkeeping work. It 
causes a total to be printed twice if the break in control is 
both minor and major. 

We conclude the discussion of the tabulator with some of 
the special wiring arrangements already mentioned in the 
text (see fig. 21). Only the unusual parts of the wiring are 
shown, and it is assumed for instance that the adding brushes 
are wired to the counters. 

Part “A” of the diagram shows the operation where 
counter (2) is to subtract algebraically numbers from the 
cards of one set from those of another. The set which is to 
be subtracted has an “x” in column 4 while both sets have 
an “x” in column 23 for negative numbers. 

Part ‘“B” of the diagram shows part of the wiring for 
adding in counter (1) only the numbers from the first card of 
each control group, but without stopping to break control 
(see ch. III, sec. 2). Columns 74-76 are used to determine 
the control groups. The wire shown on the diagram as lead- 
ing from ‘TO SHUNT PLUGS” to “PLUG TO C”’ should lead 
instead to column 17 of the ‘‘ CONTROL BRUSHES.” The cards 

should contain an “‘x”’ in column 17. 

Wiring “C”’ is for the case where special control cards are 
interspersed between the detail cards and we want one break 
in control with each control card (see ch. III, sec. 6). The 
control cards must have an ‘‘x” in column 73 and the detail 
cards not. The major control switch must be turned on and 
the minor off. 

Wiring “ D” will cause a break in control after every other 
card (see ch. V, sec. 4). For this operation the group indi- 
cation and major control switches must be on and the minor 


off. The wire labeled “to c.c.” is to be connected to the card 
count on the plugboard. 


The Reproducer 


Figure 22 shows the arrangement of the plugboard of the 
reproducer. This diagram should be compared with figure 11. 
As one would expect from figure 11, reproducing requires 
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wires from ‘‘READING BRUSHES” to “ PUNCH MAGNETS, ”’ while 
gang punching requires wires from ‘‘PUNCH BRUSHES" to 
‘“‘ PUNCH MAGNETS.” 

There are eighty ‘COMPARING RELAYS” each with two 
holes in the plugboard. In checking the work of the repro- 
ducer one of these holes must be connected to a comparing 
brush and the other to the corresponding punch brush. In 
checking the work of the gang punch one is connected to the 
reading brush and the other to a comparing brush. 

The class selectors are the same as on the tabulator. The 
two holes at the bottom are common and are for the wire to 
the ‘‘x"’ brush. 

There may be from 6 to 10 “x” brushes which precede the 
reading brushes, and from 6 to 10 more which precede the 
punch brushes. The holes in the plugboard for them are just 
above those for the punch brushes. They are used to operate 
the class selectors, and to prevent the punching of a card with 
an ‘‘x,”’ as described below. 

The ‘‘X-ELIMINATOR” is used to separate ‘x’ and “y” 
in a given column from the digits 0 to 9. For instance if we 
have x and y and the digits 0 to 9 in the same column we can 
form a new card with the x and y in one column and the 0 to 
9 in another. This can be done for five different columns at 
the same time. 

The holes marked “zero” furnish a source of “zero” 
impulses just as the card count of the tabulator gave a “1.” 

The pair of common holes “p-x”’ when wired to a ‘‘ PUNCH 
X BRUSH” will supress all punching for cards with an “x” 
(see ch. IV, sec. 3). A special switch makes it possible to 
punch only when there is or only when there is not an “x.” 


The holes ‘‘R-x” serve a similar purpose when connected to 
a ‘‘READING X BRUSH.” 


, 


The Multiplier 


The plugboard of a new type standard machine is shown in 
figure 23. This should be studied in connection with the 
description in chapter II. 
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The eighty “BRUSHES” are used to read numbers from the 
card as it is fed into the machine, These numbers may be 
read into (a) counter (1) ‘‘MULTIPLIER,” (6) counter (2) 
““sfULTIPLICAND,” (c) counter (4) and four positions of counter 
(3) “Rr. u. c.,” (d) counter (6) and four positions of counter 
(5) “L. H. c.,” and (e) counter (7) ‘‘SUMMARY COUNTER.” 

The computed products are read into counter (7) through 
the tenth row of holes ‘SUMMARY CTR.” from counters (4)-(5) 
here labeled ‘PRODUCTS COUNTER.” They are also punched 
on the card from the same counters through the twelfth row 
of holes. It is also possible to punch from counter (7) through 
the thirteenth row of holes. The wires from rows 12 and 13 
go, of course, to the eighty holes labeled ‘‘ PUNCH.” 

In order to ‘“‘round”’ a computed quantity a ‘5’’ may be 
added in any one of the sixteen positions “} ENTRY,’’ either 
for individual cards “IND.” or for each “‘GRoupP.”’ The group 
would be determined by an ‘‘x”’ punch on a card. For con- 
venience there are 6 ‘‘X-BRUSHES”’ and we can plug from any 
one of these to “c.” 

The four rows of twelve holes labeled ‘“‘ NORMAL PLUGGING” 
and “a +B—c” are associated with multiplication and 
crossfooting. For normal multiplication these holes are 
wired to connect the multiplying mechanism with the right- 
hand and left-hand components counters. Since this wiring is 
through the plugboard it is possible to change it for unusual 
operations such as forming two products at one time. 

The remainder of the holes in the plugboard are used in- 
stead of switches for controlling the machine. For instance 
if we want a fixed multiplier we wire together the two bottom 
holes in that box. If we do not want a fixed multiplier we 
must connect the top pair. Some of these switches require 
two circuits or even four. 
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Managers shall furnish a maximum estimate of cost to an 
applicant whose project is deemed acceptable. 

6. The Board of five Managers shall be appointed as 
follows: one by the International Business Machines Cor- 
poration, one by the Department of Astronomy of Columbia 
University, and three by the Advisory Council. The three 
appointed by the Advisory Council shall each serve a term of 
three years and shall not be immediately eligible for reap- 
pointment after serving two terms or portions of terms. The 
terms shall be so arranged that one of the three is appointed 
every year. A vacancy occurring before the expiration of any 
term shall be filled by the Advisory Council for the unexpired 
portion of that term. A member of the Advisory Council 
shall be eligible for membership on this Board, but if so 
appointed shall vacate his seat on the Council. 

The Board of Managers shall appoint a Chairman, Secre- 
tary, and Treasurer from its own members. The presence of 
three members shall constitute a quorum for the transaction 
of business. Majority ruling of those present shall prevail, 
the Chairman having an extra vote in the case of a tie. 

7. The President of the American Astronomical Society 
shall be requested to appoint the five members of the Advisory 
Council, each to serve for a term of five years and not to be 
eligible for immediate reappointment. The terms shall be so 
arranged that one member shall retire each year. A vacancy 
before the end of any term is to be filled by the Advisory 
Council for its unexpired portion. The Council shall elect its 
own Chairman. 

8. The duties of the Advisory Council shall be: (1) to ap- 
point three Managers as set forth in paragraph 6; (2) to give 
advice when requested by the Chairman of the Board of 
Managers; (3) to meet at least once a year with the Board 
of Managers at which meeting the Chairman of the latter shall 
report on the work of the Bureau; (4) to act as trustees of any 
funds or property acquired by the Bureau, decision as to the 
use of such property, however, to rest with the Board of 
Managers. 
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a favorable vote of six at a com- 
and Council provided notice of 
very member of both at least 


Is, 12, 116 
ic machines, % 


Ct, dhs 


Averages, 2 
Background correction, 93 
Balance control, 121 

Banks, printing, 16 

Bessel’s formula, 41. 46 
Break control, 15, 34, 49, 120 
Brushes, adding, 12, 118 
comparing, 21 
control, 13, 118 
movable, 10, 17 
reading, 4, 21, 126 
sorting, 10 

#” 18, 21 

wire, 4 


. 
: 
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Calculation control 
108 
Card count, 16, 118 
files at Bureau, 115 
gauge, 26 
punch, 8 
punching speed, 9 
Cards, for breaking 
colored, 34 
control, 34 
detail, 21, 23 
fields on, 7 
used in in 
interpretation 


Fi A i Re OG RT NE gO NRE CT a OT Nee ie ites lM eT wo nated Ok ap 


134 


Control relays, 119 

switch, calculation, 8, 22, 77, 108 
Conversion of coordinates, 81 
Correlation coefficients, 2, 31, 33 
Cos (ijh) table, 59 
Cosine series, 58, 60, 65 
Counters, 12, 18 
Counting attachment, 11 

sorter, 10 
Critical tables, 36 
Crossfeoting, 17, 19, 20, 30 


Dealing cards, 24 

Detail cards, 21, 23 

Difference formulae, 43 

Differences, 2, 39, 50, 55 

Differential equations, 75, 97, 101 

Differentiation, numerical, 43 

Digit emitter, 24 

Digiting, progressive, 31 

Disturbing function, 113 

Double punching, 10 

Duplicates, 21 

Duplicating key punch, 8 
rack, 9 

Duplicator, 8 


Eject key on card punch, 9 
Electrical multiplication table, 18 
Electromagnets, 4, 12, 116 
End figure method, 
Equations, differential, 75, 97, 101 


normal, 2, 95 


Exponential series, 65 
Feeding mechanism on tabulator, 
15 
on multiplier, 17 


yperation on reproducer, 21 


on cards, 7 


Puncuep Carp Metnop 


12, 


Gearshift lever, 121 
General perturbations, 101, 112 
Group indication, 16, 118, 122 


Half pick up, 19 
Harmonic analysis, 57, 62 
series, 2, 57, 65, 102 
High-speed listing device, 24 
reproducer, 20 
summary punch, 24 
Hollerith method, 1, 4 


Integration, numerical, 22, 43, 5 
105 
Interpolation, 41, 43 
bibliography, 43 
formulae, 41, 46 
by Lagrangian method, 47 
Interpreter, 7, 21 


Key punch, duplicating, 8, 16 
Keyboard computing machines, 25 


Lagrangian formulae, 43, 46 
method of interpolation, 47 

Least squares, 31, 33, 94 

List bank, 16 

Listing, 16 
check, 27 

Logarithms, 1, 37 

Lunary theory, 97 


Machine methods, 25 
Machines, accounting, 7, 12 
not used at Bureau, 22 


Magnitude of stars, determination, 90 


Major control, 15 
Master card, 9, 17, 21, 108 
Mechanical quadrature, 43 
verifier, 9, 26 
Minor control, 15 
Minus signs, 121 
Movable brushes, 10, 17 
Multiplicand counters, 18 
Multiplication, check, 28 
short cuts, 29 
Multiplication of series, 65, 74 
Multiplier, 17 
counters, 18 


key, duplicating. 


Punched card method, | 


Punching, double, 10 
gang, 21, 30 
off, 26 

Punching speed, 9 


ds 
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Speed, punching, 9 
reproducer, 21 
reset, 16 
sorting, 12 
summary punch, 17 
tabulator, 16 
Star catalogues, construction, 79, 92 
Stellar photometry, 90 
Stirling’s formula, 41, 46 
Subtraction tabulator, 12, 13 
Summary cards, 17, 28 
products counters, 18 
punch, 16 
punch, high-speed, 24 
Summation with tabulator, 36 
Superfluous difference method of in- 
terpolation, 54 
Switch, control, calculation, 8, 
108 
Switches, miscellaneous, 
122 


Symbol printing, 121 


on tabulator, 


Synthesis, harm: 


Tables 


subdivision of, 53 


tical 3 
critical, JO 


use in punched card form, 39 


Puncuep Carp Metruop 


Tabular functions, tables of, 36 
Tabulating card, 4, 6 

operation, 16 
Tabulation checks, 27 
Tabulator, alphabetic, 

Hollerith, 1, 12 
Test scoring machine, 
Type bars, 116 


Verifier, mechanical, 9, 26 


Watson Computing Bureau, iii, 7, 129 
Wire brushes, 4 
Wiring, 117, 124 


“x” brushes, 18, 21 
cards, 16 

control, 16 
distributor, 120 
eliminator, 126 


punch, 6, 19 


‘y” punch, 6 
Yale zone catalogue of stars, 79 


DEPARTMENT CF 


APPLIED PHYSICS 
y / b ; oe | 


UNIVERSITY OF TORONTO 


ee 


Eckert, Wallace John 
Punched card methods 
in scientific computation 


PLEASE DO NOT REMOVE 
CARDS OR SLIPS FROM THIS POCKET 


UNIVERSITY OF TORONTO LIBRARY 


